Cosmology of the Very Early Universe 



Robert H. Brandenberger 
Physics Department, McGill University, Montreal, Quebec, H3A 2T8, Canada 

In these lectures I focus on early universe models which can explain the currently observed struc- 
ture on large scales. I begin with a survey of inflationary cosmology, the current paradigm for 
understanding the origin of the universe as we observe it today. I will discuss some progress and 
problems in inflationary cosmology before moving on to a description of two alternative scenarios - 
the Matter Bounce and String Gas Cosmology. All early universe models connect to observations 
via the evolution of cosmological perturbations - a topic which will be discussed in detail in these 
lectures. 
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I. INTRODUCTION 



Overview 



Observational cosmology is currently in its "Golden 
Age". A wealth of new observational results are being 
uncovered. The cosmic microwave background (CMB) 
has been measured to high precision, the distribution of 
visible and of dark matter is being mapped out to greater 
and greater depths. New observational windows to probe 
the structure of the universe are opening up. To explain 
these observational results it is necessary to consider pro- 
cesses which happened in the very early universe. 

The inflationary scenario [J (see also 0-11]) is the cur- 
rent paradigm for the evolution of the very early uni- 
verse. Inflation can explain some of the puzzles which 
the previous paradigm - Standard Big Bang cosmology 
- could not address. More importantly, however, infla- 
tionary cosmology gave rise to the first explanation for 
the origin of inhomogeneities in the universe based on 
causal physics fsl (see also [1, 0, 0|)). The theory of 
structure formation in inflationary cosmology predicted 
the detailed shape of the angular power spectrum of CMB 
anisotropies, a prediction which was verified many years 
later observationally [1]. 

In spite of this phenomenological success, inflation- 
ary cosmology is not without its conceptual problems. 
These problems motivate the search for alternative pro- 
posals for the evolution of the early universe and for 
the generation of structure. These alternatives must be 
consistent with the current observations, and they must 
make predictions with which they can be observationally 
distinguished from inflationary cosmology. In these lec- 
tures I will discuss two alternative scenarios, the "Matter 
Bounce" (see e.g. 0) and "String Gas Cosmology" [l0| . 

I begin these lectures with a brief survey of the Stan- 
dard Big Bang (SBB) model and its problems. Possibly 
the most important problem is the absence of a structure 
formation scenario based on causal physics. I will then 
introduce inflationary cosmology and the two alternatives 
and show how they address the conceptual problems of 
the SBB model, in particular how structure is formed in 
these models. 

Since the information about the early universe is en- 
coded in the spectrum of cosmological fluctuations about 



the expanding background, it is these fluctuations which 
provide a link between the physics in the very early uni- 
verse and current cosmological observations. Section 2 is 
devoted to an overview of the theory of linearized cosmo- 
logical perturbations, the main technical tool in modern 
cosmology. The analysis of this section is applicable to all 
early universe theories. In Section 3 I give an overview of 
inflationary cosmology, focusing on the basic principles 
and emphasizing recent progress and conceptual prob- 
lems. 

The conceptual problems of inflationary cosmology 
motivate the search for alternatives. In Section 4 the 
"Matter Bounce" alternative is discussed, in Section 5 
the "String Gas" alternative. 



B. Standard Big Bang Model 

Standard Big Bang cosmology is based on three prin- 
ciples. The first is the assumption that space is homo- 
geneous and isotropic on large scales. The second is the 
assumption that the dynamics of space-time is described 
by the Einstein field equations. The third basis of the 
theory is that matter can be described as a superposi- 
tion of two classical perfect fluids, a radiation fluid with 
relativistic equation of state and pressure-less (cold) mat- 
ter. 

The first principle of Standard Cosmology implies that 
the metric takes the following most general form for a 
homogeneous and isotropic four-dimensional universe 



dt^ - a{tfdx^ 



(1) 



where t is physical time, x denote the three co-moving 
spatial coordinates (points at rest in an expanding space 
have constant co-moving coordinates) , and the scale fac- 
tor a{t) is proportional to the radiius of space. For sim- 
plicity we have assumed that the universe is spatially flat. 
The expansion rate H{t) of the universe is given by 



H{t) = - , 
a 



(2) 



where the overdot represents the derivative with respect 
to time. The cosmological redshift z{t) at time t yields 
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the amount of expansion which space has undergone be- 
tween the time t and the present time to'. 



z{t) + 1 = 



a{t) 



(3) 



In the case of a homogeneous and isotropic metric ([T|) , 
the Einstein equations which describe how matter in- 
duces curvature of space-time reduce to a set of ordinary 
differential equations, the Friedmann equations. Written 
for simplicity in the case of no cosmological constant and 
no spatial curvature they are 
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(4) 



where p is the energy density of matter and G is Newton's 
gravitational constant, and 
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P + ip) 



or equivalently 



P 



-3H{p + p). 



(5) 



(6) 



where p is the pressure density of matter. 

The third principle of Standard Cosmology says that 
the energy density and pressure are the sums of the con- 
tributions from cold matter (symbols with subscript "m" ) 
and radiation (subscripts "r" ) : 



P = Pm+ Pr 
P = Pni + Pr 



(7) 



with pm = and pr = 1/ 3pr ■ 

The first principle of Standard Cosmology, the homo- 
geneity and isotropy of space on large scales, was initially 
introduced soleley to simplify the mathematics. However, 
it is now very well confirmed observationally by the near 
isotropy of the cosmic microwave background and by the 
convergence to homogeneity in the distribution of galax- 
ies as the length scale on which the universe is probed 
increases. 

Standard Cosmology rests on three observational pil- 
lars. Firstly, it explains Hubble's redshift-distance re- 
lationship. Secondly, and most importantly for the 
theme of these lectures, it predicted the existence and 
black body nature of the cosmic microwave background 
(CMB). The argument is as follows. If we consider regu- 
lar matter and go back in time, the temperature of matter 
increases. Eventually, it exceeds the ionization temper- 
ature of atoms. Before that time, matter was a plasma, 
and space was permeated by a thermal gas of photons. A 
gas of photons will thus remain at the current time. The 
photons last scatter at a time tree, the "time of recom- 
bination" which occurs at a redshift of about 10'^. After 
that, the gas of photons remains in a thermal distribu- 
tion with a temperature T which redshifts as the universe 
expands, i.e. T ^ a^^. The CMB is this remnant gas 



FIG. 1: A sketch of the temperature (vertical axis) - time 
(horizontal axis) relation in Standard Cosmology. As the be- 
ginning of time is approached, the temperature diverges. 



of photons from the early universe. The precision mea- 
surement of the black body nature of the CMB [ll| can 
be viewed as the beginning of the "Golden Age" of ob- 
servational cosmology. The third observational pillar of 
Standard Cosmology is the good agreement between the 
predicted abundances of light elements and the observed 
ones. This agreement tells us that Standard Cosmology 
is a very good description of how the universe evolved 
back to the time of nucleosynthesis. Any modifications 
to the time evolution of Standard Cosmology must take 
place before then. 

At the present time to matter is dominated by the cold 
pressureless component. Thus, as we go back in time, 
the universe is (except close to the present time when 
it appears that a residual cosmological constant is be- 
ginning to dominate the cosmological dynamics) initially 
in a matter-dominated phase during which a{t) ^ t^/^. 
Since the energy density in cold matter scales as a"^ 
whereas that of relativistic radiation scales as a~^, there 
will be a time before which the energy density in radia- 
tion was larger than that of cold matter. The time when 
the two energy densities are equal is t^q and corresponds 
to a redshift of around z = 10^. For t < teq the uni- 
verse is radiation-dominated and a{t) ~ t^^"^. As we go 
back into the past the density and temperature increase 
without bound and a singularity is reached at which en- 
ergy density, temperature and curvature are all infinite. 
This is the "Big Bang" (see Fig. [T] for a sketch of the 
temperature/time relationship in Standard Cosmology). 



C. Problems of the Standard Big Bang 

Obviously, the assumptions on which Standard Cos- 
mology is based break down long before the singularity is 
reached. This is the "singularity problem" of the model, 
a problem which also arises in inflationary cosmology. It 
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FIG. 2: Sketch illustrating the horizon problem of Standard 
Cosmology: our past horizon at t^ec is larger than the causal 
horizon (forward light cone) at that time, thus making it im- 
possible to causally explain the observed isotropy of the CMB. 



is wrong to say that Standard Cosmology predicts a Big 
Bang. Instead, one should say that Standard Cosmology 
is incomplete and one does not understand the earliest 
moments of the evolution of the universe. 

There are, however, more "practical" "problems" of 
the Standard Big Bang model - problems in the sense 
that the model is unable to explain certain key features 
of the observed universe. The first such problem is the 
"horizon problem" [l| : within Standard Cosmology there 
is no possible explanation for the observed homogene- 
ity and isotropy of the CMB. Let us consider photons 
reaching us from opposite angles in the sky. As sketched 
in Figure m the source points for these photons on the 
last scattering surface are separated by a distance greater 
than the horizon at that time. Hence, there is no possible 
causal mechanism which can relate the temperatures at 
the two points. 

A related problem is the "size problem" : if the spatial 
sections of the universe are finite, then the only length 
scale available at the Planck time t^i is the Planck length. 
However, it we extrapolate the size of our currently ob- 
served horizon back to when the temperature was equal 
to the Planck temperature, then the corresponding wave- 
length is larger than 1/^m, many orders of magnitude 
larger than the Planck length. Standard Cosmology of- 
fers no explanation for these initial conditions. 

The third mystery of Standard Cosmology concerns 
the observed degree of spatial flatness of the universe. 
At the current time the observed energy density equals 
the "critical" energy density pc of a spatially fiat universe 
to within 10% or better. However, in Standard Cosmol- 
ogy p = Pc 'i^ an unstable critical point in an expanding 
universe. As the universe expands, the relative differ- 
ence between p and pc increases. This can be seen by 
taking the Friedmann equation in the presence of spatial 



curvature 



where 



+ eT^ = p 

3 



{aTf 



(8) 
(9) 



T being the temperature and k the curvature constant 
which is A: = ±1 or fc = for closed, open or flat spa- 
tial sections, and comparing ([8|) with the corresponding 
equation in a spatially flat universe (fc = and pc re- 
placing p. If entropy is conserved (as it is in Standard 
Cosmology) then e is constant and we obtain 



P- Pc 
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T- 



(10) 



Pc SttG Pc 

Hence, to explain the currently observed degree of spa- 
tial flatness, the initial spatial curvature had to have been 
tuned to a very high accuracy. This is the "flatness prob- 
lem". 

We observe highly non-random correlations in the dis- 
tribution of galaxies in the universe. The only force which 
can act on the relevant distances is gravity. Gravity is a 
weak force, and therefore the seed fluctuations which de- 
velop into the observed structures had to have been non- 
randomly distributed at the time teq^ the time when grav- 
itational clustering begins (see Section 2). However, as 
illustrated in Figure IH the physical length corresponding 
to the fluctuations which we observe today on the largest 
scales (they have constant comoving scale) is larger than 
the horizon at that time. Hence, there can be no causal 
generation mechanism for these perturbations ^. This is 
the "fluctuation problem" of Standard Cosmology. 

There are also more conceptual problems: Standard 
Cosmology is based on treating matter as a set of perfect 
fluids. However, we know that at high energies and tem- 
peratures a classical description of matter breaks down. 
Thus, Standard Cosmology must break down at suffi- 
ciently high energies. It does not contain the adequate 
matter physics to describe the very early universe. Sim- 
ilarly, the singularity of space-time which Standard Cos- 
mology contains corresponds to a breakdown of the as- 
sumptions on which General Relativity is based. This is 
the "singularity problem" of Standard Cosmology. 

All of the early universe scenarios which I will dis- 
cuss in the following provide solutions to the formation 
of structure problem. They can successfully explain the 
wealth of observational data on the distribution of mat- 
ter in the universe and on the CMB anisotropics. The 
rest of these lectures will focus on this point. However, 
the reader should also ask under which circumstances the 
scenarios to be discussed below address the other prob- 
lems mentioned above. 



^ This is the usual textbook argument. The student is invited to 
find (at least) two flaws in this argument. 
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FIG. 3: Sketch illustrating the formation of structure prob- 
lem of Standard Cosmology: the physical wavelength of fluc- 
tuations on fixed comoving scales which correspond to the 
large-scale structures observed in the universe is larger than 
the horizon at the time teq of equal matter and radiation, the 
time when matter fluctuations begin to grow. Hence, it is 
impossible to explain the origin of non-trivial correlations of 
the seeds for the fluctuations which had to have been present 
at that time. 
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FIG. 5: Space-time sketch of inflationary cosmology. The ver- 
tical axis is time, the horizontal axis corresponds to physical 
distance. The solid line labelled k is the physical length of 
a fixed comoving fluctuation scale. The role of the Hubble 
radius and the horizon are discussed in the text. 



t. tR to 

FIG. 4: A sketch showing the time line of inflationary cos- 
mology. The period of accelerated expansion begins at time 
ti and end at tn. The time evolution after tn corresponds to 
what happens in Standard Cosmology. 



D. Inflation as a Solution 

The idea of inflationary cosmology is to add a period 
to the evolution of the very early universe during which 
the scale factor undergoes accelerated expansion - most 
often nearly exponential growth. The time line of infla- 
tionary cosmology is sketched in Figure |4l The time ti is 
the beginning of the inflationary period, and tf; is its end 
(the meaning of the subscript R will become clear later). 
Although inflation is usually associated with physics at 
very high energy scales, e.g. E ^ lO^^Gev, all that is 
required from the initial basic considerations is that in- 
flation ends before the time of nucleosynthesis. 

During the period of inflation, the density of any pre- 
existing particles is red-shifted. Hence, if inflation is to be 
viable, it must contain a mechanism to heat the universe 
at tn, a. "reheating" mechanism - hence the subscript R 
on ti{. This mechanism must involve dramatic entropy 
generation. It is this non-adiabatic evolution which leads 
to a solution of the flatness problem, as the reader can 
verify by inspecting equation (|10l) and allowing for en- 



tropy generation at the time tn.. 

A space-time sketch of inflationary cosmology is given 
in Figure [5] . The vertical axis is time, the horizontal 
axis corresponds to physical distance. Three different 
distance scales are shown. The solid line labelled by k 
is the physical length corresponding to a flxed comoving 
perturbation. The second solid line (blue) is the Hubble 
radius 

lH{t)=H-\t). (11) 

The Hubble radius separates scales where microphysics 
dominates over gravity (sub-Hubble scales) from ones on 
which the effects of microphysics are negligible (super- 
Hubble scales) ^. Hence, a necessary requirement for a 
causal theory of structure formation is that scales we ob- 
serve today originate at sub-Hubble lengths in the early 
universe. The third length is the "horizon" , the forward 
light cone of our position at the Big Bang. The horizon 
is the causality limit. Note that because of the exponen- 
tial expansion of space during inflation, the horizon is 
exponentially larger than the Hubble radius. It is impor- 
tant not to confuse these two scales. Hubble radius and 
horizon are the same in Standard Cosmology, but in all 
three early universe scenarios which will be discussed in 
the following they are different. 



^ This statement will be demonstrated later in these lectures. 
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From Fig. [S]it is clear that provided that the period of 
inflation is sufficiently long, all scales which are currently 
observed originate as sub-B[ubble scales at the beginning 
of the inflationary phase. Thus, in inflationary cosmol- 
ogy it is possible to have a causal generation mechanism 
of fluctuations • Since matter pre-existing at ti is 
redshifted away, we are left with a matter vacuum. The 
inflationary universe scenario of structure formation is 
based on the hypothesis that all current structure orig- 
inated as quantum vacuum fluctuations. From Figure 
[5] it is also clear that the horizon problem of standard 
cosmology can be solved provided that the period of in- 
flation lasts sufficiently long. The reader should convince 
him/herself that the required period of inflation is about 
50H^^ if inflation takes place at an enegy scale of about 
lO^^GeV. Inflation thus solve both the horizon and the 
structure formation problems. 

To obtain exponential expansion of space in the con- 
text of Einstein gravity the energy density must be con- 
stant. Thus, during inflation the total energy and size 
of the universe both increase exponentially. In this way, 
inflation can solve the size and entropy problems of Stan- 
dard Cosmology. 

To summarize the main point concerning the genera- 
tion of cosmological fluctuations (the main theme of these 
lectures) in inflationary cosmology: the first crucial cri- 
terium which must be satisfied in order to have a suc- 
cessful theory of structure formation is that fiuctuation 
scales originate inside the Hubble radius. In infiationary 
cosmology it is the accelerated expansion of space dur- 
ing the inflationary phase which provides this success. In 
the following we will emphasize what is responsible for 
the corresponding success in the two alternative scenarios 
which we will discuss. 



E. Matter Bounce as a Solution 



The matter bounce scenario assumes that the contract- 
ing phase is matter-dominated at the times when scales 
we observe today exit the Hubble radius. A model in 
which the contracting phase is the time reverse of our 
current expanding phase would obey this condition. The 
assumption of an initial matter-dominated phase will be 
seen later in these lectures to be important if we want to 
obtain a scale-invariant spectrum of cosmological pertur- 
bations from initial vacuum fluctuations Il3-[l4l. 











= 1/k X 





FIG. 6: Space-time sketch in the matter bounce scenario. The 
vertical axis is conformal time rj, the horizontal axis denotes 
a co-moving space coordinate. Also, "H"^ denotes the co- 
moving Hubble radius. 



The flrst alternative to cosmological inflation as a the- 
ory of structure formation is the "matter bounce" , an 
alternative which is not yet well appreciated (for an 
overview the reader is referred to Q). The scenario is 
based on a cosmological background in which the scale 
factor a{t) bounces in a non-singular manner. 

Figure |6] shows a space-time sketch of such a bounc- 
ing cosmology. Without loss of generality we can ad- 
just the time axis such that the bounce point (minimal 
value of the scale factor) occurs at t = 0. There are 
three phases in such a non-singular bounce: the initial 
contracting phase during which the Hubble radius is de- 
creasing linearly in \t\, a bounce phase during which a 
transition from contraction to expansion takes place, and 
thirdly the usual expanding phase of Standard Cosmol- 
ogy. There is no prolonged inflationary phase after the 
bounce, nor is there a time-symmetric deflationary con- 
tracting period before the bounce point. As is obvious 
from the Figure, scales which we observe today started 
out early in the contracting phase at sub-Hubble lengths. 



Let us make a flrst comparison with inflation. A non- 
deflationary contracting phase replaces the accelerated 
expanding phase as a mechanism to bring flxed comov- 
ing scales within the Hubble radius as we go back in time, 
allowing us to consider the possibility of a causal gener- 
ation mechanism of fluctuations. Starting with vacuum 
fluctuations, a matter-dominated contracting phase is re- 
quired in order to obtain a scale- invariant spectrum. This 
corresponds to the requirement in inflationary cosmology 
that the accelerated expansion be nearly exponential. 

With Einstein gravity and matter satisfying the usual 
energy conditions it is not possible to obtain a non- 
singular bounce. However, as mentioned before, it is 
unreasonable to expect that Einstein gravity will pro- 
vide a good description of the physics near the bounce. 
There are a large number of ways to obtain a non- 
singular bouncing cosmology. To mention but a few, it 
has been shown that a bouncing cosmology results nat- 
urally from the special ghost-free higher derivative grav- 
ity Lagrangian introduced in [15]. Bounces also arise in 
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the higher-derivative non-singular universe construction 
of [l^, in "mirage cosmology" (see e.g. [l3|) which is the 
cosmology on a brane moving through a curved higher- 
dimensional bulk space-time (the time dependence for 
a brane observer is induced by the motion through the 
bulk), or - within the context of Einstein gravity - by 
making use of "quintom matter" (matter consisting of 
two components, one with regular kinetic term in the 
action, the other one with opposite sign kinetic action) 
[Tsj . For an in-depth review of ways of obtaining bounc- 
ing cosmologies see [l^ . Very recently, it has been shown 
(20.] that a bouncing cosmology can easily emerge from 
Horava-Lifshitz gravity f2l| , a new approach to quantiz- 
ing gravity. 

In the matter bounce scenario the universe begins cold 
and therefore large. Thus, the size problem of Standard 
Cosmology does not arise. As is obvious from Figure [51 
there is no horizon problem for the matter bounce sce- 
nario as long as the contracting period is long (to be 
specific, of similar duration as the post-bounce expand- 
ing phase until the present time). By the same argu- 
ment, it is possible to have a causal mechanism for gen- 
erating the primordial cosmological perturbations which 
evolve into the structures we observe today. Specifically, 
as will be discussed in the section of the matter bounce 
scenario, if the fluctuations originate as vacuum pertur- 
bations on sub-Hubble scales in the contracting phase, 
then the resulting spectrum at late times for scales exit- 
ing the Hubble radius in the matter-dominated phase of 
contraction is scale-invariant The propagation 

of infrared (IR) fluctuations through the non-singular 
bounce was analyzed in the case of the higher deriva- 
tive gravity model of fis'l in ^2^] , in mirage cosmology in 
[TtI . in the case of the quintom bounce in [H, |23| and 
for a Horava-Lifshitz bounce in 25|. The result of these 



studies is that the scale- invariance of the spectrum before 
the bounce goes persists after the bounce as long as the 
time period of the bounce phase is short compared to the 
wavelengths of the modes being considered. Note that if 
the fluctuations have a thermal origin, then the condition 
on the background cosmology to yield scale-invariance of 
the spectrum of fluctuations is different [26j . 




p = 



p - rho / 3 



FIG. 7: The dynamics of string gas cosmology. The vertical 
axis represents the scale factor of the universe, the horizon- 
tal axis is time. Along the horizontal axis, the approximate 
equation of state is also indicated. During the Hagedorn phase 
the pressure is negligible due to the cancellation between the 
positive pressure of the momentum modes and the negative 
pressure of the winding modes, after time tn the equation of 
state is that of a radiation-dominated universe. 



early quasi-static phase to the radiation phase of Stan- 
dard Cosmology. Fig. [7] shows a sketch of the evolution 
of the scale factor in string gas cosmology. 

In Figure [8] we sketch the space-time diagram in string 
gas cosmology. Since the early Hagedorn phase is quasi- 
static, the Hubble radius is infinite. For the same reason, 
the physical wavelength of fluctuations remains constant 
in this phase. At the end of the Hagedorn phase, the 
Hubble radius decreases to a microscopic value and makes 
a transition to its evolution in Standard Cosmology. 

Once again, we see that fluctuations originate on sub- 
Hubble scales. In string gas cosmology, it is the existence 
of a quasi-static phase which leads to this result. As will 
be discussed in the section on string gas cosmology, the 
source of perturbations in string gas cosmology is ther- 
mal: string therniodynamical fluctuations in a compact 
space with stable winding modes in fact leads to a scale- 
invariant spectrum (30l |. 



F. String Gas Cosmology as a Solution 



String gas cosmology [l^l (see also [11], and see [13] 
for a comprehensive review) is a toy model of cosmology 
which results from coupling a gas of fundamental strings 
to a background space-time metric. It is assumed that 
the spatial sections are compact. It is argued that the 
universe starts in a quasi-static phase during which the 
temperature of the string gas hovers at the Hagedorn 
value [2^, the maximal temperature of a gas of closed 
strings in thermal equilibrium. The string gas in this 
early phase is dominated by strings winding the compact 
spatial sections. The annihilation of winding strings will 
produce string loops and lead to a transition from the 



II. THEORY OF COSMOLOGICAL 
PERTURBATIONS 

The key tool which is used in modern cosmology to 
connect theories of the very early universe with cosmo- 
logical observations today is the theory of cosmological 
perturbations. In the following, we will give an overview 
of this theory (similar to the overview in |32l| w hich in 
turn is based on the comprehensive review [33j). We 
begin with the analysis of perturbations in Newtonian 
cosmology, a useful exercise to develop intuition and no- 
tations. 
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FIG. 8: Space-time diagram (sketch) showing the evolution of 
fixed co-moving scales in string gas cosmology. The vertical 
a:xis is time, the horizontal axis is physical distance. The solid 
curve represents the Einstein frame Hubble radius which 
shrinks abruptly to a micro-physical scale at tn and then in- 
creases linearly in time for t > tn. Fixed co-moving scales 
(the dotted lines labeled by fci and k2) which are currently 
probed in cosmological observations have wavelengths which 
are smaller than the Ifubble radius before ta. They exit the 
Hubble radius at times ti(k) just prior to tn, and propagate 
with a wavelength larger than the Hubble radius until they 
reenter the Hubble radius at times tf{k). 

A. Newtonian Theory 

The growth of density fluctuations is a consequence of 
the purely attractive nature of the gravitational force. 
Imagine (first in a non-expanding background) a density 
excess 6p localized about some point x in space. This 
fluctuation produces an attractive force which pulls the 
surrounding matter towards x. The magnitude of this 
force is proportional to Sp. Hence, by Newton's second 
law 

5p - G6p, (12) 

where G is Newton's gravitational constant. Hence, there 
is an exponential instability of flat space-time to the de- 
velopment of fluctuations. 

Obviously, in General Relativity it is inconsistent to 
consider density fluctuations in a non-expanding back- 
ground. If we consider density fluctuations in an expand- 
ing background, then the expansion of space leads to a 



friction term in (jl2p . Hence, instead of an exponential in- 
stability to the development of fluctuations, fluctuations 
in an expanding Universe will grow as a power of time. 
We will now determine what this power is and how it 
depends both on the background cosmological expansion 
rate and on the length scale of the fluctuations. 

We flrst consider the evolution of hydrodynamical mat- 
ter fluctuations in a fixed non-expanding background. In 
this context, matter is described by a perfect fluid, and 
gravity by the Newtonian gravitational potential ip. The 
fluid variables are the energy density p, the pressure p, 
the fluid velocity v, and the entropy density s. The basic 
hydrodynamical equations are 

p + V ■ (pv) = 

V + (v • V)v + -Vp + V(/3 = 
P 

S/^ip = AnGp (13) 
s + {-v ■ V)s = 

P = P{p,s). 

The flrst equation is the continuity equation, the second 
is the Euler (force) equation, the third is the Poisson 
equation of Newtonian gravity, the fourth expresses en- 
tropy conservation, and the last describes the equation 
of state of matter. The derivative with respect to time is 
denoted by an over-dot. 

The background is given by the energy density po, the 
pressure po I vanishing velocity, constant gravitational po- 
tential ipo and constant entropy density sq. As mentioned 
above, it does not satisfy the background Poisson equa- 
tion. 

The equations for cosmological perturbations are ob- 
tained by perturbiirg the fluid variables about the back- 
ground, 

p = pa + 5p 
V = 5v 

p ^ Pa + 5p (14) 
s ~ so + 5s , 

where the fluctuating fields 6p, (5v, 5p, Sip and Ss are func- 
tions of space and time, by inserting these expressions 
into the basic hydrodynamical equations (1131) , by lineariz- 
ing, and by combining the resulting equations which are 
of first order in time. We get the following differential 
equations for the energy density fluctuation 6p and the 
entropy perturbation Ss 

Sp-clv'^Sp-AnGpaSp = aV^Ss (15) 
Ss =0, 

where the variables and a describe the equation of 
state 

Sp = clSp + aSS (16) 
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with 



= (-), 



(17) 



denoting the square of the speed of sound. 

The fluctuations can be classified as foUows: If 6s 
vanishes, we have adiabatic fluctuations. If ds is non- 
vanishing but 6p = 0, we speak of an entropy fluctua- 
tion. 

The first conclusions we can draw from the basic per- 
turbation equations p3|) are that 

1) entropy fluctuations do not grow, 

2) adiabatic fluctuations are time-dependent, and 

3) entropy fluctuations seed adiabatic ones. 

All of these conclusions will remain valid in the relativis- 
tic theory. 

Since the equations are linear, we can work in Fourier 
space. Each Fourier component 6pk(t) of the fluctuation 
field Sp{x, t) evolves independently. In the case of adi- 
abatic fluctuations, the cosmological perturbations are 
described by a single field which obeys a second order 
differential equation and hence has two fundamental so- 
lutions. We will see that this conclusion remains true in 
the relativistic theory. 

Taking a closer look at the equation of motion for dp, 
we see that the third term on the left hand side represents 
the force due to gravity, a purely attractive force yield- 
ing an instability of flat space-time to the development 
of density fluctuations (as discussed earlier, see (|12l) '). 
The second term on the left hand side of (IT5|) represents 
a force due to the fluid pressure which tends to set up 
pressure waves. In the absence of entropy fluctuations, 
the evolution of 6p is governed by the combined action 
of both pressure and gravitational forces. 

Restricting our attention to adiabatic fluctuations, we 
see from that there is a critical wavelength, the Jeans 
length, whose wavenumber kj is given by 



kj = 



AttGpq.1/2 



(18) 



Fluctuations with wavelength longer than the Jeans 
length (fc <C kj) grow exponentially 



Spk{t) ^ e"^* with LUk - 4(7rGpo)'^' 



(19) 



whereas short wavelength modes (k ^ kj) oscillate with 
frequency ^ Cgk. Note that the value of the Jeans 
length depends on the equation of state of the back- 
ground. For a background dominated by relativistic ra- 
diation, the Jeans length is large (of the order of the 
Hubble radius H~^(t)), whereas for pressure-less matter 
it vanishes. 

Let us now improve on the previous analysis and study 
Newtonian cosmological fluctuations about an expand- 
ing background. In this case, the background equations 
are consistent (the non- vanishing average energy density 
leads to cosmological expansion). However, we are still 
neglecting general relativistic effects (the fluctuations of 



the metric). Such effects turn out to be dominant on 
length scales larger than the Hubble radius H^^{t), and 
thus the analysis of this section is applicable only to 
smaller scales. 

The background cosmological model is given by the en- 
ergy density po{t), the pressure Po{t), and the recessional 
velocity vq = H{t)x, where x is the Euclidean spatial 
coordinate vector ("physical coordinates"). The space- 
and time-dependent fluctuating flelds are deflned as in 
the previous section: 



p(i,x) = po{t){l + S,{t,x)) 
v{t,x) = vo(t,x) + (5v(t,x) 
p{t,x.) = po{t) +Sp{t,x.), 



(20) 



where 5^ is the fractional energy density perturbation 
(we are interested in the fractional rather than in the 
absolute energy density fluctuation!), and the pressure 
perturbation dp is defined as in (jl6p . In addition, there 
is the possibility of a non- vanishing entropy perturbation 
defined as in p^ . 

We now insert this ansatz into the basic hydrodynami- 
cal equations (jl3p . linearize in the perturbation variables, 
and combine the first order differential equations for 
and Sp into a single second order differential equation for 
dpg. The result simplifies if we work in "comoving co- 
ordinates" q which are the coordinates painted onto the 
expanding background, i.e. x(i) — a{t)q{t) . After some 
algebra, we obtain the following equation which describes 
the time evolution of density fiuctuations: 



2H6, - 



AirGpoSe 



poa^ 



(21) 



where the subscript q on the V operator indicates that 
derivatives with respect to comoving coordinates are 
used. In addition, we have the equation of entropy con- 
servation SS — . 

Comparing with the equations (jlSp obtained in the ab- 
sence of an expanding background, we see that the only 
difference is the presence of a Hubble damping term in 
the equation for Se- This term will moderate the ex- 
ponential instability of the background to long wave- 
length density fluctuations. In addition, it will lead to a 
damping of the oscillating solutions on short wavelengths. 
More specifically, for physical wavenumbers kp <^ kj 
(where kj is again given by ([H])), and in a matter- 
dominated background cosmology, the general solution 
of (|2T|) in the absence of any entropy fluctuations is given 
by 



4(t) 



Cl 



(22) 



where ci and C2 are two constants determined by the 
initial conditions, and we have dropped the subscript e 
in expressions involving S^. There are two fundamental 
solutions, the first is a growing mode with Sk{t) ^ a,{t), 
the second a decaying mode with Sk{t) ^ t^^ . On short 
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wavelengths, one obtains damped oscillatory motion: 

Sk{t) - a-^/^{t)exp{±iCsk J dt'a~\t')) . (23) 

The above analysis applies for t > t^q, when we are 
following the fluctuations of the dominant component of 
matter. In the radiation era, cold matter fluctuations 
only grow logarithmically in time since the dominant 
component of matter (the relativistic radiation) does not 
cluster and hence the term in the equation of motion for 
fluctuations which represents the gravitational attraction 
force is suppressed. In this sense, we can say that inho- 
mogeneities begin to cluster at the time t^q- 

As a simple application of the Newtonian equations 
for cosmological perturbations derived above, let us com- 
pare the predicted cosmic microwave background (CMB) 
anisotropics in a spatially flat universe with only baryonic 
matter - Model A - to the corresponding anisotropics in a 
flat Universe with mostly cold dark matter (pressure-less 
non-baryonic dark matter) - Model B. We start with the 
observationally known amplitude of the relative density 
fluctuations today (time to), and we use the fact that the 
amplitude of the CMB anisotropics on the angular scale 
9(k) corresponding to the comoving wavenumber k is set 
by the primordial value of the gravitational potential cj) - 
introduced in the following section - which in turn is re- 
lated to the primordial value of the density fluctuations 
at Hubble radius crossing (and not to its value at the 
time tree - see e.g. Chapter 17 of [33j). 

In Model A, the dominant component of the pressure- 
less matter is coupled to radiation between teq and tree, 
the time of last scattering. Thus, the Jeans length is com- 
parable to the Hubble radius. Therefore, for comoving 
galactic scales. A; ^ fcj in this time interval, and thus the 
fractional density contrast decreases as a{t)~^^^. In con- 
trast, in Model B, the dominant component of pressure- 
less matter couples only weakly to radiation, and hence 
the Jeans length is negligibly small. Thus, in Model B, 
the relative density contrast grows as a{t) between teq 
and tree- In the time interval tree < t < to, the fluctua- 
tions scale identically in Models A and B. Summarizing, 
we conclude, working backwards in time from a fixed 
amplitude of Sk today, that the amplitudes of Sk{teq) in 
Models A and B (and thus their primordial values) are 
related by 

Skiteq)\A - {^%^f'Skiteq)\B. (24) 

Hence, in Model A (without non-baryonic dark matter) 
the CMB anisotropics are predicted to be a factor of 
about 30 larger [sS] than in Model B, way in excess of the 
recent observational results. This is one of the strongest 
arguments for the existence of non-baryonic dark matter 
^. Note that the precise value of the enhancement factor 



In my opinion no proposed alternative to particle dark matter 



depends on the value of the cosmological constant A - the 
above result holds for A = 0. 



B. Observables 

Let us consider perturbations on a fixed comoving 
length scale given by a comoving wavenumber k. The 
corresponding physical length increases as a{t). This is 
to be compared with the Hubble radius H^-^{t) which 
scales as t provided a{t) grows as a power of t. In the late 
time Universe, a{t) ^ t^^^ in the radiation-dominated 
phase (i.e. for t < teq), and a{t) ^ t^^^ in the matter- 
dominated period (teq < t < to)- Thus, we see that at 
sufficiently early times, all comoving scales had a physi- 
cal length larger than the Hubble radius. If we consider 
large cosmological scales (e.g. those corresponding to the 
observed CMB anisotropics or to galaxy clusters), the 
time tff(fc) of "Hubble radius crossing" (when the phys- 
ical length was equal to the Hubble radius) was in fact 
later than teq- 

Cosmological fluctuations can be described either in 
position space or in momentum space. In position 
space, we compute the root mean square mass fluctu- 
ation 6M/M{k,t) in a sphere of radius I = 27r/fc at time 
t. A scale-invariant spectrum of fluctuations is defined 
by the relation 

^ (fc, iff (fc)) - const.. (25) 

Such a spectrum was first suggested by Harrison [35j and 
Zeldovich j36| as a reasonable choice for the spectrum of 
cosmological fluctuations. We can introduce the "spec- 
tral index" n of cosmological fluctuations by the relation 

{^)\k,tHm ^ k-\ (26) 

and thus a scale- invariant spectrum corresponds to n = 1. 

To make the transition to the (more frequently used) 
momentum space representation, we Fourier decompose 
the fractional spatial density contrast 

(5e(x,t) = V^/^ j d^kSeik,t)e'^ '' , (27) 

where V is a. cutoff volume. The povirer spectrum Pg 

of density fluctuations is deflned by 

Ps{k) ^ k^\Se{k)\\ (28) 

where k is the magnitude of k, and we have assumed 
for simplicity a Gaussian distribution of fluctuations in 



should be taken serious unless their proponents can demonstrate 
that their model can reproduce the agreement between the ampli- 
tudes of the CMB anisotropies on one hand and of the large-scale 
matter power spectrum on the other hand. 
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which the amplitude of the fluctuations only depends on 
k. 

We can also define the power spectrum of the gravita- 
tional potential (p: 

P^(fc) - fc3|^~^(fc)|2. (29) 

These two power spectra are related by the Poisson equa- 
tion (fT3| 

P^{k) ^ k-^Ps{k). (30) 

In general, the condition of scale-invariance is ex- 
pressed in momentum space in terms of the power spec- 
trum evaluated at a fixed time. To obtain this condition, 
we first use the time dependence of the fractional density 
fluctuation from (|22p to determine the mass fluctuations 
at a fixed time t > tnik) > teq (the last inequality is a 
condition on the scales considered) 

The time of Hubble radius crossing is given by 

a{tH{k))k-^ = 2tH{k), (32) 

and thus 

tHikf'^^k-^. (33) 

Inserting this result into (|3T|) and making use of (|26|) we 
find 

(^)^(fc,t) ^ (34) 

Since, for reasonable values of the index of the power 
spectrum, dM/M{k,t) is dominated by the Fourier 
modes with wavenumber fc, we find that (j34]) implies 

\Se\^-^k", (35) 

or, equivalently, 

P^(fc)~fc"-i. (36) 

C. Classical Relativistic Theory 

1. Introduction 

The Newtonian theory of cosmological fluctuations dis- 
cussed in the previous section breaks down on scales 
larger than the Hubble radius because it neglects per- 
turbations of the metric, and because on large scales the 
metric fluctuations dominate the dynamics. 

Let us begin with a heuristic argument to show why 
metric fluctuations are important on scales larger than 
the Hubble radius. For such inhomogeneities, one should 
be able to approximately describe the evolution of the 



space-time by applying the first FRW equation (|571) of 
homogeneous and isotropic cosmology to the local Uni- 
verse (this approximation is made more rigorous in |37^ ) . 
Based on this equation, a large-scale fluctuation of the en- 
ergy density will lead to a fluctuation ( "(5a" ) of the scale 
factor a which grows in time. This is due to the fact that 
self gravity amplifies fluctuations even on length scales A 
greater than the Hubble radius. 

This argument is made rigorous in the following anal- 
ysis of cosmological fluctuations in the context of general 
relativity, where both metric and matter inhomogeneities 
are taken into account. We will consider fluctuations 
about a homogeneous and isotropic background cosmol- 
ogy, given by the metric ([1}, which can be written in 
conformal time i] (defined by dt = a(t)dri) as 

ds^ = air]f{dr,^ -dx") . (37) 

The theory of cosmological perturbations is based on 
expanding the Einstein equations to linear order about 
the background metric. The theory was initially devel- 
oped in pioneering works by Lifshitz [ssj . Significant 
progress in the understanding of the physics of cosmo- 
logical fluctuations was achieved by Bardeen [s^ who re- 
alized the importance of subtracting gauge artifacts (see 
below) from the analysis (see also [40,|41[). The following 
discussion is based on Part I of the comprehensive review 
article j33| . Other reviews - in some cases emphasizing 
different approaches - are [4^ - |45| . 

2. Classifying Fluctuations 

The flrst step in the analysis of metric fluctuations is to 
classify them according to their transformation proper- 
ties under spatial rotations. There are scalar, vector and 
second rank tensor fluctuations. In linear theory, there is 
no coupling between the different fluctuation modes, and 
hence they evolve independently (for some subtleties in 
this classification, see j46|). 

We begin by expanding the metric about the FRW 
background metric gljS given by (|37p: 

5m- = 9l°J + • (38) 

The background metric depends only on time, whereas 
the metric fluctuations Sg^^, depend on both space and 
time. Since the metric is a symmetric tensor, there are 
at first sight 10 fluctuating degrees of freedom in ^g^^. 

There are four degrees of freedom which correspond 
to scalar metric fluctuations (the only four ways of con- 
structing a metric from scalar functions): 

= a' ( 2(^PS^E,,) ) ' ^39) 

where the four fluctuating degrees of freedom are de- 
noted (following the notation of [s^) (j),B,E, and ip, a 
comma denotes the ordinary partial derivative (if we had 
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included spatial curvature of the background metric, it 
would have been the covariant derivative with respect to 
the spatial metric), and 6ij is the Kronecker symbol. 

There are also four vector degrees of freedom of metric 
fluctuations, consisting of the four ways of constructing 
metric fluctuations from three vectors: 

^■'^^ - ( A, ) ' (40) 

where Si and Fi are two divergence- less vectors (for a vec- 
tor with non- vanishing divergence, the divergence con- 
tributes to the scalar gravitational fluctuation modes). 

Finally, there are two tensor modes which correspond 
to the two polarization states of gravitational waves: 

= -a^ ( ) , (41) 
where hij is trace-free and divergence-less 

K^Hj^O. (42) 

Gravitational waves do not couple at linear order to 
the matter fluctuations. Vector fluctuations decay in an 
expanding background cosmology and hence are not usu- 
ally cosmologically important. The most important fluc- 
tuations, at least in inflationary cosmology, are the scalar 
metric fluctuations, the fluctuations which couple to mat- 
ter inhomogeneities and which are the relativistic gener- 
alization of the Newtonian perturbations considered in 
the previous section. Note that vector and tensor per- 
turbations are important (see e.g. [l^) in topological 
defects models of structure formation and in bouncing 
cosmologies (48l |. 

3. Gauge Transformation 

The theory of cosmological perturbations is at flrst 
sight complicated by the issue of gauge invariance (at 
the flnal stage, however, we will see that we can make 
use of the gauge freedom to substantially simplify the 
theory). The coordinates t,x of space-time carry no in- 
dependent physical meaning. They are just labels to des- 
ignate points in the space-time manifold. By performing 
a small-amplitude transformation of the space-time co- 
ordinates (called "gauge transformation" in the follow- 
ing), we can easily introduce "flctitious" fluctuations in 
a homogeneous and isotropic Universe. These modes are 
"gauge artifacts" . 

We will in the following take an "active" view of gauge 
transformation. Let us consider two space-time man- 
ifolds, one of them a homogeneous and isotropic Uni- 
verse A4q, the other a physical Universe Ai with inho- 
mogeneities. A choice of coordinates can be considered 
to be a mapping D between the manifolds Mo and M. 
Let us consider a second mapping V which will take the 
same point (e.g. the origin of a flxed coordinate system) 



in A^o to a different point in A4. Using the inverse of 
these maps V and V, we can assign two different sets of 
coordinates to points in A4. 

Consider now a physical quantity Q (e.g. the Ricci 
scalar) on M, and the corresponding physical quantity 
(5*-°^ on TMo Then, in the flrst coordinate system given 
by the mapping V, the perturbation 6Q of Q at the point 
p E Ai is deflned by 

6Q{p) ^ Q{p) - Q("\V-\p)) . (43) 

Analogously, in the second coordinate system given by 
I?, the perturbation is deflned by 

5g(p) = Q(p)-g(")(p-i(p)). (44) 

The difference 

AQ(p) = S~Qip) ~ 5Q{p) (45) 

is obviously a gauge artifact and carries no physical sig- 
niflcance. 

Some of the degrees of freedom of metric perturbation 
introduced in the flrst subsection are gauge artifacts. To 
isolate these, we must study how coordinate transfor- 
mations act on the metric. There are four independent 
gauge degrees of freedom corresponding to the coordinate 
transformation 

^ £^ = a;*' + . (46) 

The zero (time) component of leads to a scalar 
metric fluctuation. The spatial three vector can be 
decomposed as 

e = a+7% (47) 

(where 7*-' is the spatial background metric) into a trans- 
verse piece Q.^ which has two degrees of freedom which 
yield vector perturbations, and the second term (given by 
the gradient of a scalar f ) which gives a scalar fluctuation. 
To summarize this paragraph, there are two scalar gauge 
modes given by and ^, and two vector modes given by 
the transverse three vector Q.^. Thus, there remain two 
physical scalar and two vector fluctuation modes. The 
gravitational waves are gauge- invariant. 

Let us now focus on how the scalar gauge transforma- 
tions (i.e. the transformations given by and ^) act 
on the scalar metric fluctuation variables (f>, E, and ip. 
An immediate calculation yields: 

a 

B = B + i°-^' (48) 
E = E-^ 

^ - ^ + 

a 

where a prime indicates the derivative with respect to 
conformal time rj. 
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There are two approaches to deal with the gauge ambi- 
guities. The first is to fix a gauge, i.e. to pick conditions 
on the coordinates which completely eliminate the gauge 
freedom, the second is to work with a basis of gauge- 
invariant variables. If one wants to adopt the gauge-fixed 
approach, there are many different gauge choices. Note 
that the often used synchronous gauge determined by 
ggOfi _ Q (Joes not totally fix the gauge. A convenient 
system which completely fixes the coordinates is the so- 
called longitudinal or conformal Newtonian gauge 
defined by B ^ E = Q. This is the gauge we will use in 
the following. 

Note that the above analysis was strictly at linear order 
in perturbation theory. Beyond linear order, the struc- 
ture of perturbation theory becomes much more involved. 
In fact, one can show [4^ that the only fluctuation vari- 
ables which are invariant under all coordinate transfor- 
mations are perturbations of variables which are constant 
in the background space-time. 

4. Equations of Motion 

We begin with the Einstein equations 

G^. = SttGT^. , (49) 

where G^j/ is the Einstein tensor associated with the 
space-time metric g^^, and is the energy- momentum 
tensor of matter, insert the ansatz for metric and matter 

perturbed about a FRW background (5^S'('7), V^'^Hi]))- 

5^^(x,77) = 9^^}{il) + 5g^,^{yL,iq) (50) 
(p(x,77) = </Jo(?7) + (5(p(x,77) , (51) 

(where we have for simplicity replaced general matter by 
a scalar matter field ^p) and expand to linear order in the 
fluctuating fields, obtaining the following equations: 

<5G^. = S^GJT^.. (52) 

In the above, Sg^^u is the perturbation in the metric and 
5(f is the fluctuation of the matter field (f. 
Inserting the longitudinal gauge metric 

ds^ = [(1 + 2(j))dTf - (1 - 2iP)-i,jdx''dx^] (53) 

into the general perturbation equations ([5^ yields the 
following set of equations of motion: 

- 3H(-H(/> + V') + V^^/; = A-KGa^ST^ 

[{2H' +H'^)(f> + H4'' +^" +2H'4j']6} (54) 
+ ^V^D5]~^f''D.kj = -AnGa^Sr;, 

where D = (j) — and % — a' / a. 



The first conclusion we can draw is that if no 
anisotropic stress is present in the matter at linear or- 
der in fiuctuating fields, i.e. ^Tj = for i 7^ j, then the 
two metric fluctuation variables coincide: 

(j) = ^. (55) 

This will be the case in most simple cosmological models, 
e.g. in theories with matter described by a set of scalar 
fields with canonical form of the action, and in the case 
of a perfect fluid with no anisotropic stress. 

Let us now restrict our attention to the case of matter 
described in terms of a single scalar field f (in which 
case the perturbations on large scales are automatically 
adiabatic) which can be expanded as 

V3(x,77) = Lpa{ii) + 5p{yi,ri) (56) 

in terms of background matter Lp^ and matter fiuctuation 
(593(x, 77). Then, in longitudinal gauge, ([51]) reduce to the 
following set of equations of motion 



(I)' +H(f> = 
<l)" +3n(f)' + {H +2H'^)(f, = 



4:nG{(p'oSip' + Va'^Sp) 
AnGip'^Sip (57) 
A7rG{ip'QSip' - v'a^Sip) , 



where V denotes the derivative of V with respect to if. 
These equations can be combined to give the following 
second order differential equation for the relativistic po- 
tential (j): 

0" + 2 1 -H - ^ 1 - v^<j) + 2\n - 10 = 0. 

(58) 

This is the final result for the classical evolution of 
cosmological fiuctuations. First of all, we note the sim- 
ilarities with the equation (PT|) obtained in the Newto- 
nian theory. The final term in ([55)1 is the force due to 
gravity leading to the instability, the second to last term 
is the pressure force leading to oscillations (relativistic 
since we are considering matter to be a relativistic field) , 
and the second term represents the Hubble friction. For 
each wavenumber there are two fundamental solutions. 
On small scales (fc > H), the solutions correspond to 
damped oscillations, on large scales (fc < H) the oscil- 
lations freeze out and the dynamics is governed by the 
gravitational force competing with the Hubble friction 
term. Note, in particular, how the Hubble radius natu- 
rally emerges as the scale where the nature of the fluctu- 
ating modes changes from oscillatory to frozen. 

Considering the equation in a bit more detail, we ob- 
serve that if the equation of state of the background is 
independent of time (which will be the case if "H — — 
0), then in an expanding background, the dominant mode 
of (|58)) is constant, and the sub-dominant mode decays. 
If the equation of state is not constant, then the dominant 
mode is not constant in time. Specifically, at the end of 
inflation H < 0, and this leads to a growth of cj) (see the 
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following subsection). In contrast, in a contracting phase 
the dominant mode of (f> is growing. 

To study the quantitative impUcations of the equation 
of motion (|58|) . it is convenient to introduce jSOt l5l| the 
variable C (which, up to a correction term of the order 
V^^ which is unimportant for large-scale fluctuations, is 
equal to the curvature perturbation TZ in comoving gauge 
[5l)by 
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where 



1 



P 

w = — 
P 



(59) 



(60) 



characterizes the equation of state of matter. In terms of 
C, the equation of motion ((55)) takes on the form 



-CHil + w) 



0(V^ 



(61) 



On large scales, the right hand side of the equation is 
negligible, which leads to the conclusion that large-scale 
cosmological fluctuations satisfy 



C{l+w) = 0. 



(62) 



This implies that C is constant except possibly if 1+w = 
at some point in time during the cosmological evolution 
(which occurs during reheating in inflationary cosmology 
if the inflaton field undergoes oscillations [53] )■ In sin- 
gle matter field models it is indeed possible to show that 
C = on super- Hubble scales independent of assumptions 
on the equation of state [13, ■ This "conservation law" 
makes it easy to relate initial fluctuations to final fluctu- 
ations in inflationary cosmology, as will be illustrated in 
the following subsection. 

In the presence of non-adiabatic fluctuations (entropy 
fluctuations) there is a source term on the righ-hand side 
of (|62p which is proportional to the amplitude of the en- 
tropy fluctuations. Thus, as already seen in the case of 
the Newtonian theory of cosmological fluctuations, non- 
adiabatic fluctuations induce a growing adiabatic mode 
(see [H, [131 for discussions of the consequences in double 
field inflationary models). 



D. Quantum Theory 

1. Overview 

As we have seen in the first section, in many models 
of the very early Universe, in particular in inflationary 
cosmology, in string gas cosmology and in the matter 
bounce scenario, primordial inhomogeneities are gener- 
ated in an initial phase on sub-Hubble scales. The wave- 
length is then stretched relative to the Hubble radius, 
becomes larger than the Hubble radius at some time and 



then propagates on super-Hubble scales until re-entering 
at late cosmological times. In a majority of the cur- 
rent structure formation scenarios (string gas cosmology 
is an exception in this respect), fluctuations are assumed 
to emerge as quantum vacuum perturbations. Hence, 
to describe the generation and evolution of the inhomo- 
geneities a quantum treatment is required. 

In the context of a Universe with a de Sitter phase, 
the quantum origin of cosmological fluctuations was first 
discussed in and also [1, 0] for earlier ideas. In par- 
ticular, Mukhanov 's'l and Press realized that in an 
exponentially expanding background, the curvature fluc- 
tuations would be scale- invariant, and Mukhanov pro- 
vided a quantitative calculation which also yielded the 
logarithmic deviation from exact scale-invariance. 

The role of the Hubble radius has already been men- 
tioned repeatedly in these lectures. In particular, in the 
previous subsection we saw that the Hubble radius sep- 
arates scales on which fluctuations oscillate (sub-Hubble 
scales) from those where they are frozen in (super-Hubble 
scales) . Another way to see the role of the Hubble radius 
is to consider the equation of a free scalar matter held (p 
on an unperturbed expanding background: 



V2 

</3 + 3Hip -ip 



0. 



(63) 



The second term on the left hand side of this equation 
leads to damping of ip with a characteristic decay rate 
given hy H. As a consequence, in the absence of the 
spatial gradient term, ip would be of the order of mag- 
nitude Hip. Thus, comparing the second and the third 
terms on the left hand side, we immediately see that the 
microscopic (spatial gradient) term dominates on length 
scales smaller than the Hubble radius, leading to oscil- 
latory motion, whereas this term is negligible on scales 
larger than the Hubble radius, and the evolution of ip is 
determined primarily by gravity. Note that in general 
cosmological models the Hubble radius is much smaller 
than the horizon (the forward light cone calculated from 
the initial time). In an inflationary universe, the hori- 
zon is larger by a factor of at least exp(A^), where N 
is the number of e-foldings of inflation, and the lower 
bound is taken on if the Hubble radius and horizon coin- 
cide until inflation begins. It is very important to realize 
this difference, a difference which is obscured in most ar- 
ticles on cosmology in which the term "horizon" is used 
when "Hubble radius" is meant. Note, in particular, that 
the homogeneous inflaton fleld contains causal informa- 
tion on super-Hubble but sub-horizon scales. Hence, it is 
completely consistent with causality [11] to have a micro- 
physical process related to the background scalar matter 
fleld lead to exponential ampliflcation of the amplitude 
of fluctuations during reheating on such scales, as it does 
in models in which entropy perturbations are present and 
not suppressed during inflation [H, 113] ■ Note that also in 
string gas cosmology and in the matter bounce scenario 
the Hubble radius and horizon are completely different. 
There are general relativistic conservation laws (ssj 
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which imply that adiabatic fluctuations produced locally 
must be Poisson-statistic suppressed on scales larger than 
the Hubble radius. For example, fluctuations produced 
by the formation of topological defects at a phase tran- 
sition in the early universe are initially isocurvature (en- 
tropy) in nature (see e.g. for a discussion). Via the 
source term in the equation of motion (jisp , a growing adi- 
abatic mode is induced, but at any fixed time the spec- 
trum of the curvature fluctuation on scales larger than 
the Hubble radius has index n = 4 (Poisson). A similar 
conclusion applies to models [HO, El of modulated re- 
heating (see J32J for a nice discussion), and to models in 
which moduli fields obtain masses after some symmetry 
breaking, their quantum fluctuations then inducing cos- 
mological perturbations. A prototypical example is given 
by axion fluctuations in an inflationary universe (see e.g. 
[63| and references therein). 

To understand the generation and evolution of fluctu- 
ations in current models of the very early Universe, we 
need both Quantum Mechanics and General Relativity, 
i.e. quantum gravity. At first sight, we are thus faced 
with an intractable problem, since the theory of quantum 
gravity is not yet established. We are saved by the fact 
that today on large cosmological scales the fractional am- 
plitude of the fluctuations is smaller than 1. Since gravity 
is a purely attractive force, the fluctuations had to have 
been - at least in the context of an eternally expanding 
background cosmology - very small in the early Universe. 
Thus, a linearized analysis of the fluctuations (about a 
classical cosmological background) is self-consistent. 

From the classical theory of cosmological perturbations 
discussed in the previous section, we know that the analy- 
sis of scalar metric inhomogeneities can be reduced - after 
extracting gauge artifacts - to the study of the evolution 
of a single fluctuating variable. Thus, we conclude that 
the quantum theory of cosmological perturbations must 
be reducible to the quantum theory of a single free scalar 
field which we will denote by v. Since the background in 
which this scalar field evolves is time-dependent, the mass 
off will be time-dependent. The time-dependence of the 
mass will lead to quantum particle production over time 
if we start the evolution in the vacuum state for v. As we 
will see, this quantum particle production corresponds to 
the development and growth of the cosmological fluctu- 
ations. Thus, the quantum theory of cosmological fluc- 
tuations provides a consistent framework to study both 
the generation and the evolution of metric perturbations. 
The following analysis is based on Part II of [33l] . 



2. Outline of the Analysis 

In order to obtain the action for linearized cosmological 
perturbations, we expand the action to quadratic order 
in the fluctuating degrees of freedom. The linear terms 
cancel because the background is taken to satisfy the 
background equations of motion. 

We begin with the Einstein-Hilbert action for gravity 



and the action of a scalar matter held (for the more com- 
plicated case of general hydrodynamical fluctuations the 
reader is referred to ,33]) 



S = 



J d^x^[-^R + - Vi^)] , (64) 



where R is the Ricci curvature scalar. 

The simplest way to proceed is to work in longitudinal 
gauge, in which the metric and matter take the form 
(assuming no anisotropic stress) 

ds^ = a2(?7)[(l-|-2(/)(77,x))d?72 - (l-20(i,x))dx2] 
93(77, x) = ifiaiv) + S'fiiv,^) ■ (65) 

The two fluctuation variables (f) and tp must be linked by 
the Einstein constraint equations since there cannot be 
matter fluctuations without induced metric fluctuations. 

The two nontrivial tasks of the lengthy [s^ compu- 
tation of the quadratic piece of the action is to find out 
what combination of ip and gives the variable v in terms 
of which the action has canonical kinetic term, and what 
the form of the time-dependent mass is. This calculation 
involves inserting the ansatz (I55|) into the action (|64l) . 
expanding the result to second order in the fiuctuating 
fields, making use of the background and of the constraint 
equations, and dropping total derivative terms from the 
action. In the context of scalar field matter, the quan- 
tum theory of cosmological fiuctuations was developed 
by Mukhanov [64, 65] and Sasaki [6^. The result is the 
following contribution S^^^ to the action quadratic in the 
perturbations: 



1 J d^x[^ 



v.iV,i H V 

z 



(66) 



where the canonical variable v (the "Sasaki-Mukhanov 
variable" introduced in [65| - see also ^3]) given by 



with H = a' /a, and where 



n 



(67) 



(68) 



As long as the equation of state does not change over 
time) 



z{t]) - a{ri) 



(69) 



Note that the variable v is related to the curvature per- 
turbation TZ in comoving coordinates introduced in [s^] 
and closely related to the variable ( used in [E^, HH : 



V = zTZ . 



(70) 



The equation of motion which follows from the action 
(|66p is (in momentum space) 



+k Vk- —Vk = , 



(71) 
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where Vk is the k'th Fourier mode of v. The mass term 
in the above equation is in general given by the Hubble 
scale. Thus, it immediately follows that on small length 
scales, i.e. for k > kn, the solutions for Vk are constant 
amplitude oscillations . These oscillations freeze out at 
Hubble radius crossing, i.e. when k = kn- On longer 
scales {k kn), there is a mode of Vk which scales as z. 
This mode is the dominant one in an expanding universe, 
but not in a contracting one. 

Given the action (|66p. the quantization of the cos- 
mological perturbations can be performed by canonical 
quantization (in the same way that a scalar matter field 
on a fixed cosmological background is quantized [gsI]). 

The final step in the quantum theory of cosmological 
perturbations is to specify an initial state. Since in in- 
flationary cosmology all pre-existing classical fluctuations 
are red-shifted by the accelerated expansion of space, one 
usually assumes (we will return to a criticism of this point 
when discussing the trans-Planckian problem of inflation- 
ary cosmology) that the field v starts out at the initial 
time ti mode by mode in its vacuum state. Two ques- 
tions immediately emerge: what is the initial time 
and which of the many possible vacuum states should be 
chosen. It is usually assumed that since the fluctuations 
only oscillate on sub-Hubble scales, the choice of the ini- 
tial time is not important, as long as it is earlier than the 
time when scales of cosmological interest today cross the 
Hubble radius during the inflationary phase. The state 
is usually taken to be the Bunch-Davies vacuum (see e.g. 
[68]), since this state is empty of particles at ti in the 
coordinate frame determined by the FRW coordinates 
Thus, we choose the initial conditions 



1 



(72) 



2k 

, \fk 
Vkiv.) - ^ 

where rji is the conformal time corresponding to the phys- 
ical time ti. 

Returning to the case of an expanding universe, the 
scaling 



Vk 



(73) 



implies that the wave function of the quantum variable 
Vk which performs quantum vacuum fluctuations on sub- 
Hubble scales, stops oscillating on super-Hubble scales 
and instead is squeezed (the amplitude increases in con- 
flguration space but decreases in momentum space). This 
squeezing corresponds to quantum particle production. 
It is also one of the two conditions which are required 
for the classicalization of the fluctuations. The second 
condition is decoherence which is induced by the non- 
linearities in the dynamical system which are inevitable 
since the Einstein action leads to highly nonlinear equati- 
ions (see [6^ for a recent discussion of this point, and [70j| 
for related work). 



Note that the squeezing of cosmological fluctuations 
on super-Hubble scales occurs in all models, in particu- 
lar in string gas cosmology and in the bouncing universe 
scenario since also in these scenarios perturbations prop- 
agate on super-Hubble scales for a long period of time. In 
a contracting phase, the dominant mode of Vk on super- 
Hubble scales is not the one given in (I75|) (which in this 
case is a decaying mode), but the second mode which 
scales as z~p with an exponent p which is positive and 
whose exact value depends on the background equation 
of state. 

Applications of this theory in inflationary cosmology, 
in the matter bounce scenario and in string gas cosmol- 
ogy will be considered in the respective sections of these 
lecture notes. 



3. Quantum Theory of Gravitational Waves 

The quantization of gravitational waves parallels the 
quantization of scalar metric fluctuations, but is more 
simple because there are no gauge ambiguities. Note 
that at the level of linear fluctuations, scalar metric fluc- 
tuations and gravitational waves are independent. Both 
can be quantized on the same cosmological background 
determined by the background scale factor and the back- 
ground matter. However, in contrast to the case of scalar 
metric fluctuations, the tensor modes are also present in 
pure gravity (i.e. in the absence of matter). 

Starting point is the action (|64p . Into this action we 
insert the metric which corresponds to a classical cosmo- 
logical background plus tensor metric fluctuations: 



ds^ = 0^(77) [dry^ - {5ij + hij)dx^dx^] 



(74) 



where the second rank tensor hij(r],x.) represents the 
gravitational waves, and in turn can be decomposed as 



(75) 



into the two polarization states. Here, e^- and efj are 
two fixed polarization tensors, and /i+ and are the 
two coefficient functions. 

To quadratic order in the fluctuating fields, the action 
consists of separate terms involving and h^- Each 
term is of the form 



leading to the equation of motion 



i—hk + k^hk = 0. 



(76) 



(77) 



The variable in terms of which the action (|76p has canon- 
ical kinetic term is 



^.k 



ahh 



(78) 
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and its equation of motion is 

+ = 0. (79) 

This equation is very similar to the corresponding equa- 
tion dTTj) for scalar gravitational inhomogeneities, except 
that in the mass term the scale factor a{ri) replaces z{ri), 
which leads to a very different evolution of scalar and 
tensor modes during the reheating phase in inflationary 
cosmology during which the equation of state of the back- 
ground matter changes dramatically. 

Based on the above discussion we have the following 
theory for the generation and evolution of gravitational 
waves in an accelerating Universe (first developed by Gr- 
ishchuk (7l|): waves exist as quantum vacuum fluctu- 
ations at the initial time on all scales. They oscillate 
until the length scale crosses the Hubble radius. At that 
point, the oscillations freeze out and the quantum state 
of gravitational waves begins to be squeezed in the sense 
that 

Mfc(^)^a(^), (80) 

which, from (|78p corresponds to constant amplitude of 
hk ■ The squeezing of the vacuum state leads to the emer- 
gence of classical properties of this state, as in the case 
of scalar metric fluctuations. 



III. INFLATIONARY COSMOLOGY 
A. Mechanism of Inflation 

The idea of inflationary cosmology is to assume that 
there was a period in the very early Universe during 
which the scale factor was accelerating, i.e. a > 0. This 
implies that the Hubble radius was shrinking in comoving 
coordinates, or, equivalently, that fixed comoving scales 
were "exiting" the Hubble radius. In the simplest mod- 
els of inflation, the scale factor increases nearly expo- 
nentially. As illustrated in Figure [HI the basic geome- 
try of inflationary cosmology provides a solution of the 
fluctuation problem. As long as the phase of inflation 
is sufficiently long, all length scales within our present 
Hubble radius today originate at the beginning of infla- 
tion with a wavelength smaller than the Hubble radius 
at that time. Thus, it is possible to create perturbations 
locally using physics obeying the laws of special relativity 
(in particular causality). As will be discussed later, it is 
quantum vacuum fluctuations of matter fields and their 
associated curvature perturbations which are responsible 
for the structure we observe today. 

Postulating a phase of inflation in the very early uni- 
verse solves the horizon problem of the SBB, namely it 
explains why the causal horizon at the time tree when 
photons last scatter is larger than the radius of the past 
light cone at tree, the part of the last scattering surface 
which is visible today in CMB experiments. Inflation also 




FIG. 9; Space-time diagram (sketch) showing the evolution 
of scales in inflationary cosmology. The vertical axis is time, 
and the period of inflation lasts between ti and ta, and is 
followed by the radiation-dominated phase of standard big 
bang cosmology. During exponential inflation, the Hubble 
radius is constant in physical spatial coordinates (the 

horizontal axis), whereas it increases linearly in time after tR. 
The physical length corresponding to a fixed comoving length 
scale labelled by its wavenumber k increases exponentially 
during inflation but increases less fast than the Hubble radius 
(namely as t^^^), after inflation. 



explains the near flatness of the universe: in a decelerat- 
ing universe spatial flatness is an unstable fixed point of 
the dynamics, whereas in an accelerating universe it be- 
comes an attractor. Another important aspect of the in- 
flationary solution of the flatness problem is that inflation 
exponentially increases the volume of space. Without in- 
flation, it is not possible that a Planck scale universe at 
the Planck time evolves into a sufficiently large universe 
today. 



B. Fluctuations in Inflationary Cosmology 

We will now use the quantum theory of cosmological 
perturbations developed in the previous section to calcu- 
late the spectrum of curvature fluctuations in inflationary 
cosmology. The starting point are quantum vacuum ini- 
tial conditions for the canonical fluctuation variable Vk'. 



VkiVi) = 




for all k for which the wavelength is smaller than the 
Hubble radius at the initial time ti. 

The amplitude remains unchanged until the modes exit 



17 



the Hubble radius at the respective times tnik) given by 

a-\tH{k))k = H. (82) 

We need to compute the power spectrum Vuik) of the 
curvature fluctuation TZ defined in ([70]) at some late time 
t when the modes are super-Hubble. We first relate the 
power spectrum via the growth rate (j73p of v on super- 
Hubble scales to the power spectrum at the time tnik) 
and then use the constancy of the amplitude of v on sub- 
Hubble scales to relate it to the initial conditions (fSTI) . 
Thus 



k^z-^{t)\v,{t)\^ 

'a{tH[k)) 
ez-^{tH{k))\vk{tH{k))\'' 



(83) 
VkitHm? 



a-\tH{k))\vk{Uf 



where in the final step we have used (IMl) and the con- 
stancy of the amplitude of v on sub-Hubble scales. Mak- 
ing use of the condition for Hubble radius crossing, 
and of the initial conditions (I8ip . we immediately see that 



(84) 



and that thus a scale invariant power spectrum with am- 
plitude proportional to results, in agreement with 
what was argued on heuristic grounds in the overview of 
inflation in the the first section. To obtain the precise 
amplitude, we need to make use of the relation between 
z and a. We obtain 



Vnik,t) ^ — 



(85) 



which for any given value of k is to be evaluated at the 
time tnik) (before the end of infiation). For a scalar field 
potential (see following subsection) 



(86) 



the resulting amplitude in ([85|) is A. Thus, in order to 
obtain the observed value of the power spectrum of the 
order of 10~^°, the coupling constant A must be tuned to 
a very small value. 



determined by the Friedmann-Robertson- Walker (FRW) 
equations 



a 



SttGp 



(87) 



and 



a 
a 



A-kG 



where for simplicity we have omitted the contributions 
of spatial curvature (since spatial curvature is diluted 
during infiation) and of the cosmological constant (since 
any small cosmological constant which might be present 
today has no effect in the early Universe since the asso- 
ciated energy density does not increase when going into 
the past). From (|88|) it is clear that in order to obtain an 
accelerating universe, matter with sufficiently negative 
pressure 



P<-\p 



(89) 



is required. Exponential inffation is obtained for p — —p. 

Conventional perfect ffuids have positive semi-definite 
pressure and thus cannot yield infiation. Quantum field 
theory can come to the rescue. We know that a descrip- 
tion of matter in terms of classical perfect ffuids must 
break down at early times. An improved description of 
matter will be given in terms of quantum fields. Matter 
which we observe today consists of spin 1/2 and spin 1 
fields. Such fields cannot yield inffation in the context 
of the renormalizable quantum field theory framework. 
The existence of scalar matter fields is postulated to ex- 
plain the generation of mass via spontaneous symmetry 
breaking. Scalar matter fields (denoted by (f) are spe- 
cial in that they allow - at the level of a renormalizable 
action - the presence of a potential energy term V{lp). 
The energy density and pressure of a scalar ffeld Lp with 
canonically normalized action ^ 



C 



(90) 



(where Greek indices are space-time indices and g is the 
determinant of the metric) are given by 



C. Models of Inflation 

Let us now consider how it is possible to obtain a phase 
of cosmological infiation. We will assume that space-time 
is described using the equations of General Relativity 
^. In this case, the dynamics of the scale factor a{t) is 



P 
P 



V{p) 
V{p). 



(91) 



Thus, it is possible to obtain an almost exponentially ex- 
panding universe provided the scalar ffeld configuration 



Note, however, that the first model of exponential expansion of 
space 0] made use of a higher derivative gravitational action. 



See ]72 
terms. 



for discussions of fields with non-canonical kinetic 
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^ satisfies 

i(Vp^)2 « ^(^), (92) 

« V{^). (93) 

In the above, Vp = a^^V is the gradient with respect 
to physical as opposed to comoving coordinates. Since 
spatial gradients redshift as the universe expands, the 
first condition will (for single scalar field models) always 
be satisfied if it is satisfied at the initial time ^. It is the 
second condition which is harder to satisfy. In particular, 
this condition is in general not preserved in time even it 
is initially satisfied. 

In his original model [l[, Guth assumed that infiation 
was generated by a scalar field p> sitting in a false vacuum 
with positive vacuum energy. Hence, the conditions (j92p 
and (l93l) are automatically satisfied. Inflation ends when 
Lp tunnels to its true vacuum state with (by assumption) 
vanishing potential energy. However, as Guth realized 
himself, this model is plagued by a "graceful exit" prob- 
lem: the nucleation produces a bubble of true vacuum 
whose initial size is microscopic j?!] (see also [t^ for a re- 
view of field theory methods which are useful in infiation- 
ary cosmology). After inflation, this bubble will expand 
at the speed of light, but it is too small to explain the 
observed size of the universe. Since space between the 
bubbles expands exponentially, the probability of bub- 
ble percolation is negligible, at least in Einstein gravity. 
For this reason, the focus in inflationary model building 
shifted to the "slow-roll" paradigm [TTj]. 

It is sufficient to obtain a period of cosmological in- 
fiation that the slow-roll conditions for (p are satisfied. 
Recall that the equation of motion for a homogeneous 
scalar field in a cosmological space-time is (as follows 
from dini)) is 

>p + 3Hip = -V'ip), (94) 

where a prime indicates the derivative with respect to ip. 
In order that the scalar field roll slowly, it is necessary 
that 

<p < 3Hip (95) 

such that the first term in the scalar field equation of 
motion ((94)) is negligible. In this case, the condition ((93)) 
becomes 

iy)^ < 487rG (96) 

and ((^ becomes 

V" 

— < 247rG. (97) 



^ The scalar field yielding inflation is called the inflaton. 
In fact, careful studies show that since the gradients de- 
crease even in a non-inflationary backgrounds, they can become 
subdominant even if they arc not initially subdominant. 



There are many models of scalar field-driven slow-roll 
inflation. Many of them can be divided into three groups: 
small-field infiation, large-field inflation and hybrid infla- 
tion. Small-field inflationary models are based on ideas 
from spontaneous symmetry breaking in particle physics. 
We take the scalar field to have a potential of the form 

V{v) = Ixi^p' - <j'r : (98) 

where a can be interpreted as a symmetry breaking scale, 
and A is a dimensionless coupling constant. The hope of 
initial small-field models ("new infiation" was that 
the scalar field would begin rolling close to its symmetric 
point p — 0, where thermal equilibrium initial conditions 
would localize it in the early universe. At sufficiently high 
temperatures, = is a stable ground state of the one- 
loop finite temperature effective potential Vriip) (see e.g. 
(Tq for a review) . Once the temperature drops to a value 
smaller than the critical temperature Tc, p> — turns into 
an unstable local maximum of Vt{p), and p is free to roll 
towards a ground state of the zero temperature potential 
((M)) . The direction of the initial rolling is triggered by 
quantum fluctuations. The reader can easily check that 
for the potential ((98p the slow-roll conditions cannot be 
satisfied if cr rupi , where m^; is the Planck mass which 
is related to G. If the potential is modified to a Coleman- 
Weinberg [tII form 

F(^) = V[lnM_i]+i^A.^ (99) 
4 V 4 16 

(where v denotes the value of the minimum of the poten- 
tial) then the slow-roll conditions can be satisfied. How- 
ever, this corresponds to a severe fine-tuning of the shape 
of the potential. A further problem for most small-field 
models of infiation (see e.g. 79] for a review) is that in 
order to end up close to the slow-roll trajectory, the initial 
field velocity must be constrained to be very small. This 
initial condition problem of small-field models of inflation 
effects a number of recently proposed brane infiation sce- 
narios, see e.g. [s^ for a discussion. 

There is another reason for abandoning small-field in- 
fiation models: in order to obtain a sufficiently small 
amplitude of density fiuctuations, the interaction coeffi- 
cients of (p must be very small. This makes it inconsistent 
to assume that (p started out in thermal equilibrium [sj . 
In the absence of thermal equilibrium, the phase space 
of initial conditions is much larger for large values of ip. 

This brings us to the discussion of large-field inflation 
models, initially proposed in [82l| under the name "chaotic 
inflation" . The simplest example is provided by a massive 
scalar field with potential 

Vi^p) = imV' , (100) 

where m is the mass. It is assumed that the scalar field 
rolls towards the origin from large values of |(/?|. It is a 
simple exercise for the reader to verify that the slow-roll 
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conditions (|M|) and ([97)) are satisfied provided 

1 



z-ITLpl 



(101) 



Values oi\(p\ comparable or greater than rupi are also re- 
quired in other realizations of large- field inflation. Hence, 
one may worry whether such a toy model can consis- 
tently be embedded in a realistic particle physics model, 
e.g. supergravity. In many such models V{(p) receives 
supergravity-induced correction terms which destroy the 
flatness of the potential for \(p\ > uipi. As can be seen 
by applying the formulas of the previous subsection to 
the potential pOOp . a value of m ^ lO^^GeV is required 
in order to obtain the observed amplitude of density 
fluctuations. Hence, the configuration space range with 
\(p\ > nipi but V{lp) < m^i dominates the measure of 
field values. It can also be verified that the slow-roll tra- 
jectory is a local attractor in field initial condition space 
[74| . even including metric fiuctuations at the perturba- 
tive level jssj . 

With two scalar fields it is possible to construct a class 
of models which combine some of the nice features of 
large-field infiation (large phase space of initial condi- 
tions yielding inflation) and of small-field inflation (bet- 
ter contact with conventional particle physics). These 
are models of hybrid inflation 84] . To give a prototyp- 
ical example, consider two scalar fields tp and x with a 
potential 



1 



1 



IgW- (102) 



V{ip,x) = -Ax(x^-a^)^ + -mV 

In the absence of thermal equilibrium, it is natural to 
assume that \ip\ begins at large values, values for which 
the effective mass of x is positive and hence x begins 
at X = 0. The parameters in the potential (|102[) are 
now chosen such that if is slowly rolling for values of 
\ip\ somewhat smaller than Wpj, but that the potential 
energy for these field values is dominated by the first 
term on the right-hand side of (|102p . The reader can 
easily verify that for this model it is no longer required 
to have values of |(^| greater than rUpi in order to obtain 
slow-rolling ^ The field (p is slowly rolling whereas the 
potential energy is determined by the contribution from 
X- Once \ip\ drops to the value 



9 



(103) 



the configuration x = becomes unstable and decays to 
its ground state |x| = a, yielding a graceful exit from 
infiation. Since in this example the ground state of x is 
not unique, there is the possibility of the formation of 
topological defects at the end of infiation (see [ssl - ISTj for 
reviews of topological defects in cosmology). 



* Note that the slow-roll conditions 1 196 I I and I I97II were derived 
assuming that H is given by the contribution of ip to V which is 
not the case here. 



D. Reheating in Inflationary Cosmology 

1. Preheating 

After the slow-roll conditions break down, the period 
of infiation ends, and the inflaton ip begins to oscillate 
around its ground state. Due to couplings of ip to other 
matter flelds, the energy of the universe, which at the end 
of the period of inflation is stored completely in (p, gets 
transferred to the matter flelds of the particle physics 
Standard Model. Initially, the energy transfer was de- 
scribed perturbatively ^88j. Later, it was realized (89l - [9^ 
that through a parametric resonance instability, parti- 
cles are very rapidly produced, leading to a fast energy 
transfer ("preheating"). The quanta later thermalize, 
and thereafter the universe evolves as described by SBB 
cosmology. 

Let us give a brief overview of the theory of reheating 
in inflationary cosmology. For an in-depth review the 
student is referred to [93']. We assume that the infla- 
ton ip is coupled to another scalar field x take the 
interaction Lagrangian to be 



r 1 2 2 2 

i-int = f X : 



(104) 



where g is a dimensionless coupling constant. 

In the initial perturbative analysis of reheating (ssj 
the decay of the coherently oscillating inflaton field cp 
is treated perturbatively. The interaction rate F of pro- 
cesses in which two ip quanta interact to produce a pair of 
X particles is taken as the decay rate of the infiaton. The 
infiaton dynamics is modeled via an effective equation 



ip + 3Hip + Tip 



-Viip). 



(105) 



For small coupling constant, the interaction rate F is typ- 
ically much smaller than the Hubble parameter at the end 
of infiation. Thus, at the beginning of the phase of infia- 
ton oscillations, the energy loss into particles is initially 
negligible compared to the energy loss due to the expan- 
sion of space. It is only once the Hubble expansion rate 
decreases to a value comparable to F that x particle pro- 
duction becomes effective. It is the energy density at the 
time when H = T which determines how much energy 
ends up in x particles and thus determines the "reheat- 
ing temperature" , the temperature of the SM fields after 
energy transfer. 



Tr ^ (Tnipi) 



1/2 



(106) 



Since F is proportional to the square of the coupling con- 
stant g which is generally very small, perturbative reheat- 
ing is slow and produces a reheating temperature which 
can be very low compared to the energy scale at which 
infiation takes place. 

There are two main problems with the perturbative de- 
cay analysis described above. First of all, even if the infia- 
ton decay were perturbative, it is not justified to use the 
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heuristic equation (|105|) since it violates the fluctuation- 
dissipation theorem: in systems with dissipation, there 
are always fluctuations, and these are missing in (11051) . 
For an improved effective equation of motion see e.g. [94| . 

The main problem with the perturbative analysis is 
that it does not take into account the coherent nature of 
the inflaton field ^p. At the beginning of the period of 
oscillations </? is not a superposition of free asymptotic 
single inflaton states, but rather a coherently oscillating 
homogeneous field. The large amplitude of oscillation 
implies that it is well justified to treat the inflaton clas- 
sically. However, the matter fields can be assumed to 
start off in their vacuum state (the red-shifting during 
the period of inflation will remove any matter particles 
present at the beginning of inflation). Thus, matter fields 
X must be treated quantum mechanically. The improved 
approach to reheating initiated in [s^ (see also j95|) is 
to consider reheating as a quantum production of x par- 
ticles in a classical Lp background. 

We expand the field x is terms of the usual creation 
and annihilation operators. The mode functions satisfy 
the equation 



Xk 



g'^{tYsin'{mt)\xk = 



(107) 

where $(<:) is the amplitude of oscillation of ip. 

The first level of approximation [89] is to neglect the 
expansion of space, to estimate the efficiency of reheating 
in this approximation, and to check self-consistency of 
the approximation. The equation (jl07p then reduces to 

Xk + + ml + g2$2 sin2 [mt)) Xfe - , (108) 

and $(t ) is constant. The equation is the Mathieu equa- 
tion 19611 which is conventionally written in the form 



Xfc + (Afc-2gcos2z)xfc 



0, 



(109) 



where we have introduced the dimensionless time variable 
z = mt and a prime now denotes the derivative with 
respect to z. Comparing the coefficients, we see that 



Ak 



2q 
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The equation (I109P describes parametric resonance in 
classical mechanics. From the examples in classical me- 
chanics in which this equation arises we know that there 
are bands of k values (depending on the frequency m of 
the external source) which exhibit exponential instabil- 
ity. For the general theory of the Mathieu equation the 
reader is referred to [9^. For values of fc in a resonance 
band, the growth of Xfc can be written as 



Xk 
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where fik is called the Floquet exponent. In the model we 
are considering in this subsection, resonance occurs for all 
long wavelength modes - roughly speaking all modes with 



This is called broad-band resonance [90|,|92|. The Floquet 
exponent is of order 1, and this implies very efficient en- 
ergy transfer from the coherently oscillating infiaton field 
to a gas of X particles. This initial stage of energy transfer 
is called preheating. The time scale of the energy transfer 
is short compared to the cosmological time scale at 
the end of infiation [89j . Hence, the approximation of 
neglecting the expansion of space is self-consistent. 

It is possible to perform an improved analysis which 
keeps track of the expansion of space. The first step 
is to rescale the field variable to extract the cosmological 
red- shifting 



Xk{t) = a{t)xk{t).. 



(113) 



and to work in terms of conformal time. The field Xk 
obeys an equation without cosmological damping term 
which is very similar to (|109l) , except that the oscillatory 
correction term to the mass is replaced by a more general 
periodic function of the conformal time rj. The Floquet 
theory of [96| also applies to the resulting equation, and 
we reach the conclusion that the exponential instability 
of X (modulated by a{t)) persists when including the ef- 
fects of the expansion of space. For details the reader is 
referred to 1921 or to the recent review [9311. 



2. Preheating of Metric Perturbations 

A parametric resonance instability should be expected 
for all fields which couple to the oscillating inflaton con- 
densate. In particular, metric fluctuations should also 
be effected. As was shown in [53i] . there is no paramet- 
ric resonance instability for long wavelength modes in 
the case of purely adiabatic perturbations (see also ^3f\ ) . 
However, in the presence of entropy fluctuations, para- 
metric amplification of cosmological perturbations driven 
by the oscillating inflaton condensate is possible [H, [s^l 
(see also 

To study the preheating of metric fluctuations, we re- 
fer back to the formalism of cosmological perturbations 
developed in Subsection (|II C 4p . In the presence of an en- 
tropy perturbation SS, the "conservation" equation (j62p 
for the variable C (which describes the curvature pertur- 
bation in comoving gauge) gets replaced by 



c 



-SS. 



P- 
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k < grn$ . 



(112) 



If long wavelength fluctuations of a matter field x which 
coupled to the infiaton ip (x corresponds to an entropy 
mode) undergo exponential growth during reheating as 
they do in the case of a model with broad parametric 
resonance, then this will induce an exponential growth of 
SS which will via (|114l) induce an exponentially growing 
curvature perturbation. 

For an elegant formalism to compute the source term 
in (11141) in terms of the scalar fields (p and x the reader is 
referred to ^981 . This resonant growth of entropy fluctu- 
ations is only important in models in which the entropy 
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fluctuations are not suppressed du ring inflation. 



I during inflation. Some 
recent examples were studied in [H, llOd ] . The source for 
this instabihty need not be the osciflations of the inflaton 
fleld. In SUSY models, th e dec ay of flat directions can 
also induce this instability lOlj . 



of a hierarchy in scales, namely [1& 



(115) 



3. Effects of Noise on Reheating 

The author of these lecture notes at some point became 
worried that noise effects might eliminate the paramet- 
ric resonance instability during reheating. One source of 
noise is the quantum fluctuations in the inflaton Lp which 
might lead to non-periodic corrections to the fluctuating 
mass term which drives the resonance. However, under 
specific assumptions on the properties of the noise (the 
most important of these being that there are no correla- 
tions in the noise between different oscillation periods of 
the background) it was shown [102, 103] that the presence 
of the noise increas es the effi ciency of the resonance. 

The analysis of |l02l . Il03j | is based on mapping the 
question into a problem in random matrix theory. The 
matrices involved are the transfer matrices in phase space 
corresponding to the dynamical system. Making use of 
a theorem in random matrix theory called the "Fursten- 
berg Theorem" it can be shown that for each mode fc, the 
Floquet exponent ^j. in the presence of noise is strictly 
larger than in the absence of noise. This implies that 
the stability bands for the dynamical system disappear. 
Mapping this classical dynamical systems problem into a 
time-independent one-dimensional non-relativistic quan- 
tum mechanics problem via the standard correspondence 
yields a ne w pr oof of Anderson localization in one space 
dimension |104| . 



E. Problems of Inflation 

In spite of the phenomenological success of the infla- 
tionary paradigm, conventional scalar field-driven infla- 
tion suffers from several important conceptual problems. 

The first problem concern the nature of the infiaton, 
the scalar field which generates the inflationary expan- 
sion. No particle corresponding to the excitation of a 
scalar fleld has yet been observed in nature, and the 
Higgs fleld which is introduced to give elementary par- 
ticles masses in the Standard Model of particle physics 
does not have the required flatness of the potential to 
yield inflation, unless it is non-minimally coupled to grav- 
ity [105]. In particle physics theories beyond the Stan- 
dard Model there are often many scalar flelds, but it is 
in general very hard to obtain the required flatness prop- 
erties on the potential 

The second problem (the amplitude problem) re- 
lates to the amplitude of the spectrum of cosniological 
perturbations. In a wide class of inflationary models, ob- 
taining the correct amplitude requires the introduction 



where A(/j is the change in the inflaton field during the 
minimal length of the inflationary period, and V{if) is 
the potential energy during inflation. 

A more ser ious problem is the trans-Planckian 
problem jl07l |. Returning to the space-time diagram 
of inflation (see Figure fTO|). we can immediately deduce 
that, provided that the period of inflation lasted suffi- 
ciently long (for GUT scale inflation the number is about 
70 e- foldings), then all scales inside the Hubble radius to- 
day started out with a physical wavelength smaller than 
the Planck scale at the beginning of inflation. Now, the 
theory of cosmological perturbations is based on Ein- 
stein's theory of General Relativity coupled to a simple 
semi-classical description of matter. It is clear that these 
building blocks of the theory are inapplicable on scales 
comparable and smaller than the Planck scale. Thus, the 
key successful prediction of inflation (the theory of the 
origin of fluctuations) is based on suspect calculations 
since new physics must enter into a correct computation 
of the spectrum of cosmological perturbations. The key 
question is as to whether the predictions obtained using 
the current theory are sensitive to the specifics of the 
unknown theory which takes over on small scales. 

One approach to study the sensitivity of the usual 
predictions of infiationary cosmology to the unknown 
physics on trans-Planckian scales is to study toy models 
of ultraviolet physics which allow e xplic i t cal culations. 
The first approach which was used [lOSl Il09j | is to re- 
place the usual linear dispersion relation for the Fourier 
modes of the fluctuations by a modifled dispersion re- 
lation, a dispersion relation which is linear for physical 
wavenumbers smaller than the scale of new physics, but 
deviates on larger scales. Such dispersion relations were 
used previously to test the sensitivity of black hole radi- 
ation on th e un known physics of the UV [llOl llll| . It 
was found |l08l | that if the evolution of modes on the 
trans-Planckian scales is non-adiabatic, then substantial 
deviations of the spectrum of fluctuations from the usual 
results are possible. Non-adiabatic evolution turns an 
initial state minimizing the energy density into a state 
which is excited once the wavelength becomes larger than 
the cutoff scale. Back-reaction effects of these excitations 
may limit the magnitude of the trans- Planckian effects if 
we assume that inflation took place 112l4ll4 |. On the 
other hand, large trans-Planckian effects may also pre- 
vent the onset of inflation. 

A fourth problem is the singularity problem. It was 
known for a long time that standard Big Bang cosmol- 
ogy cannot be the complete story of the early universe 
because of the initial singularity, a singularity which is 
unavoidable when basing cosmology on Einstein's field 
equations in the presence of a m atter source obeying the 
weak energy conditions (see e.g. [llSj for a textbook dis- 
cussion). Recently, the singularity theorems have been 
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FIG. 10: Space-time diagram (sketch) of inflationary cos- 
mology including the two zones of ignorance - sub-Planckian 
wavelengths and trans-Planckian densities. The symbols have 
the same meaning as in Figure 1. Note, specifically, that - 
as long as the period of inflation lasts a couple of e-foldings 
longer than the minimal value required for inflation to address 
the problems of standard big bang cosmology - all wavelengths 
of cosmological interest to us today start out at the beginning 
of the period of inflation with a wavelength which is in the 
zone of ignorance. 



early universe. 

A final problem which we will mention here is the con- 
cern that the energy scale at which inflation takes place is 
too high to justify an effective field theory analysis based 
on Einstein gravity. In simple toy models of infiation, 
iuShe energy scale during the period of infiation is about 
lO^^GeV, very close to the string scale in many string 
models, and not too far from the Planck scale. Thus, 
correction terms in the effective action for matter and 
gravity may already be important at the energy scale of 
inflation, and the cosmological dynamics may be rather 
different from what is obtained when neglecting the cor- 
rection terms. 

In Figure [TT] we show once again the space-time sketch 
of inflationary cosmology. In addition to the length scales 
which appear in the previous versions of this flgure, we 
have now shaded the "zones of ignorance" , zones where 
the Einstein gravity effective action is sure to break down. 
As described above, fluctuations emerge from the short 
distance zone of ignorance (except if the period of infla- 
tion is very short), and the energy scale of inflation might 
put the period of inflation too close to the high energy 
density zone of ignorance to trust the predictions based 
on using the Einstein action. 

The arguments in this subsection provide a motivation 
for considering alternative scenarios of early universe cos- 
mology. Below we will focus on two scenarios, the matter 
bounce and string gas cosmology. It is important to em- 
phasize that these are not the only alternat ive s cenarios. 
Others include the Pre-B ig-Bang scenario [ll8| and the 
Ekpyrotic paradigm jllOl ]. 



generalized to apply to Einstein gravity coupled to scalar 
fleld ma tter, i.e. to scalar fleld-driven inflationary cos- 
mology |116|. It was shown that, in this context, a past 
singularity at some point in space is unavoidable. Thus 
we know, from the outset, that scalar field-driven infla- 
tion cannot be the ultimate theory of the very early uni- 
verse. 

The Achilles heel of scalar field-driven inflationary cos- 
mology may be the cosmological constant problem. 

We know from observations that the large quantum vac- 
uum energy of fleld theories does not gravitate today. 
However, to obtain a period of inflation one is using the 
part of the energy-momentum tensor of the scalar fleld 
which looks like the vacuum energy. In the absence of 
a solution of the cosmological constant problem it is un- 
clear whether scalar fleld-driven inflation is robust, i.e. 
whether the mechanism which renders the quantum vac- 
uum energy gravitationally inert today will not also pre- 
vent the vacuum energy from gravitating during the pe- 
riod of slow-rolling of the inflaton fleld. Note that the 
approach to addressing the cosmological constant prob- 
lem making use of the gravitational back-reaction of long 
range fluctuations (see [ll7j | for a summary of this ap- 
proach) does not prevent a long period of inflation in the 



IV. MATTER BOUNCE 

The matter bounce scenario is a non-singular cosmol- 
ogy in which time runs from — cxd to +oo. Negative times 
correspond to a contracting phase, positive times to ex- 
pansion. It is assumed that the bounce phase is short 

9 

As indicated in Figure |6l the fluctuations which we 
observe today have exited the Hubble radius early in the 
contracting phase. The word matter in matter bounce 
should emphasize that we are assuming that the phase 
during which the fluctuations which we observe today are 
exiting the Hubble radius is dominated by cold pressure- 
less matter. Viewing the contracting phase as the time 
reverse of the expanding phase we are living in, this last 
assumption we have made does not seem very restrictive. 
In fact, one could expect that some entropy is generated 
during the bounce, in which case the matter phase during 



^ Short means short compared to the time it takes light to travel 
over distances of about 1mm, the physical length which corre- 
sponds to to the largest comoving scales we observe today (as- 
suming that the energy scale of the bounce is about lO^^GeV). 
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FIG. 11: Space-time diagram (sketch) of inflationary cos- 
mology including the two zones of ignorance - sub-Planckian 
wavelengths and trans-Planckian densities. The symbols have 
the same meaning as in Figure 1. Note, specifically, that - 
as long as the period of inflation lasts a couple of e-foldings 
longer than the minimal value required for inflation to address 
the problems of standard big bang cosmology - all wavelengths 
of cosmological interest to us today start out at the beginning 
of the period of inflation with a wavelength which is in the 
zone of ignorance. 



the contracting period would last up to higher energy 
densities than it does in the expanding phase. 



A. Models for a Matter Bounce and Background 
Cosmology 

It follows from t he si ngularity theorems of General Rel- 
ativity (see e.g. jll5| ) that in order to obtain a non- 
singular bounce we must either invoke matter which vi- 
olates the weak energy condition, or go beyond General 
Relativity. A review of ways to obtain a non-singular 
cosmology is given in (Toj . In the following we will only 
mention a few recent attempts which the author of these 
lecture notes has been involve d in. 

Introducing quintom matter |l2Clj | is a way of obtaining 
a non-singular bouncing cosmology, as discussed in [l8| . 
Quintom matter is a set of two matter fields, one of them 
regular matter (obeying the weak energy condition) , the 
second a field with opposite sign kinetic ternij^ a field 
which violates the energy conditions. We can [18] model 



both matter components with scalar fields. Let us denote 
the mass of the regular one {(p) by m, and by M that of 
the field with wrong sign kinetic term. We assume that 
early in the contracting phase both fields are oscillating, 
but that the amplitude A of (p greatly exceeds the corre- 
sponding amplitude A of ip such that the energy density 
is dominated by p. Both fields will initially be oscillating 
during the contracting phase, and both amplitudes grow 
at the same rate. At some point, A will become so large 
that the oscillations of (p freeze out Then, A will grow 
only slowly, whereas A will continue to increase. Thus, 
the (negative) energy density in p will grow in absolute 
values relative to that of p. The total energy density will 
decrease towards 0. At that point, H = by the Fried- 
mann equations. It can in fact easily be seen that H > 
when H = 0. Hence, a non-singular bou nce occurs. The 
Higgs sector of the Lee- Wick model |l2l| provides a con- 
crete realization of the quintom bounce model, as studied 
in . Quintom models like all other models with nega- 
tive sign kinetic terms suffer from an instability problem 
[l22j | in the matter sector and are hence problematic. 

Turning to the second way of obtaining a non-singular 
bounce, namely by modifying the gravitational action, 
we should emphasize that modifications of the gravita- 
tional action are expected at the high densities at which 
the bounce will occur. Concrete examples were already 
mentioned in the Introduction: the higher derivative La- 
grangian resulting from the nonsingular universe con- 
struction of JJi], the model of which is based on 
a non-local higher derivative action which is ghost -free 
about Minkowski space-time, mirage cosmologies |l23l | 
resulting from the effective action of gravity on a brane 
which is moving into and out of a high-curvature throat 
in a higher-dimensional space-time. 

More recently, interest has focused on a non-singular 
bouncing cosmology which emerges from Horava-Lifshitz 
gravity [2l| . As shown in [23| , if the spatial sections have 
a non-vanishing spatial curvature constant fc, then the 
higher-derivative terms in the action will lead to terms 
in the Friedmann equations which act as ghost radiation 
and ghost anisotropic stress, i.e. terms of gravitational 
origin and negative effective energy density which scale 
as and a~^, respectively. Starting with a contract- 
ing universe dominated by regular matter, eventually the 
ghost terms will catch up and yield a non-singular bounce 
in analogy to how the ghost matter in the quintom model 
does. 

The analysis of the spectrum of cosmological perturba- 
tions on scales of current observational interest is, how- 
ever, independent of the details of the bouncing phase, 
as long as that phase is short compared to the time it 
takes light to travel over the length scales of current in- 
terest. We now turn to a discussion of the evolution of 
fluctuations in a generic matter bounce model. 



This corresponds to entering a region of large-field inflation. 
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Fluctuations in a Matter Bounce 



The equation of motion for the Fourier mode Vk of v is 



vl + {k^ - —)vk 



0. 



(116) 



If the equation of state of the background is time- 
independent, then z ^ a and hence the negative square 
mass term in (I116P is H^. Thus, on length scales smaller 
than the Hubble radius, the solutions of (I116P are oscil- 
lating, whereas on larger scales they are frozen in, and 
their amplitude depends on the time evolution of z. 

In the case of an expanding universe the dominant 
mode of v scales as z. However, in a contracting universe 
it is the second of the two modes which dominates. If the 
contracting phase is matter-dominated, i.e. a{t) ~ t^^^ 
and r]{t) ~ t^^^ the dominant mode of v scales as r]~^ 
and hence 



cirf' +C2r] ^ 



(117) 



where ci and C2 are constants. The ci mode is the mode 
for which C is constant on super-Hubble scales. However, 
in a contracting universe it is the C2 mode which domi- 
nates and leads to a scale- invariant spectrum [T2| - [l^ : 

Pc(fc,ry) ~ k^\vkiv)\'a-\rj) (118) 

rj ' 

^ const , 

where we have made use of the scaling of the domi- 
nant mode of Wfe, the Hubble radius crossing condition 
ilnik) ^ k~^ , and the assumption that we have a vac- 
uum spectrum at Hubble radius crossing. 

Up to this point, the analysis shows that in the 
contracting phase the curvature fluctuations are scale- 
invariant. In non-singular bouncing cosmologies, the fluc- 
tuations can be followed in an unambiguous way through 
the bounce. This was done in the case of the higher 
derivative bounce model of [l^ in [l^], and in a non- 
singular bouncing mirage cosmology model the analysis 
was performed in [l7[. In the case of the Quintom and 
Lee- Wick bounces, respectively, the fluctuations were fol- 
lowed through the bounce in [23| and respectively. 
Finally, in the case of the HL bounce, the evolution of 
the fluctuations through the bounce was recently studied 
in [25^. 

The equations of motion can be solved numerically 
without approximation. Alternatively, we can solve them 
approximately using analytical techniques. Key to the 
analytical analysis are the General Relativistic matching 
conditions for fl uctuation s across a phase transition in 
the background [l24 Il25j | . These conditions imply that 
both <I> and C, are conserved at the bounce, where 



However, as stressed in [l26i |. these matching conditions 
can only be used at a transition when the background 
metric obeys the matching conditions. This is not the 
case if we were to match directly between the contracting 
matter phase and the expanding matter phase. 

In the case of a non-singular bounce we have three 
phases: the initial contracting phase with a fixed equa- 
tion of state (e.g. u; = 0), a bounce phase during which 
the universe smoothly transits between contraction and 
expansion, and finally the expanding phase with constant 
w. We need to apply the matching conditions twice: first 
at the transition between the contracting phase and the 
bounce phase (on both sides of the matching surface the 
universe is contracting), and then between the bounc- 
ing phase and the expanding phase. The bottom line of 
the studies of [13, 12a - l25| is that on length scales large 
compared to the time of the bounce, the spectrum of 
curvature fiuctuations is not changed during the bounce 
phase. Since typically the bounce time is set by a mi- 
crophysical scale whereas the wavelength of fluctuations 
which we observe today is macroscopic (about 1mm if the 
bounce scale is set by the particle physics GUT scale), 
we conclude that for scales relevant to current observa- 
tions the spectrum is unchanged during the bounce. This 
completes the demonstration that a non-singular matter 
bounce leads to a scale-invariant spectrum of cosmologi- 
cal perturbations after the bounce provided that the ini- 
tial spectrum on sub-Hubble scales is vacuum. Initial 
thermal fiuctuations were followed through the bounce 
in [2^. Note that perturbations are processed during 
a bounce and this implies that even if the background 
cosmology were cyclic, the full evolution including linear 
fluctuations would be non-cyclic 31 1. 



C = C+T7 



1 p$ 
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C. Key Predictions of the Matter Bounce 

Ganonical single field inflation models predict very 
small non-Gaussianities in the spectrum of fluctuations. 
One way to characterize the non-Gaussianities is via the 
three point function of the curvature flu ctua tion, also 
called the "bispectrum" . As realized in (l27{ . the bis- 
pectrum induced in the matter bounce scenario has an 
amplitude which is at the borderline of what the Planck 
satellite experiment will be able to detect, and it has a 
special form. These are specific predictions of the matter 
bounce scenario with which the matter bounce scenario 
can be distinguished from those of standard infiationary 
models (see [128] for a recent detailed review of non- 
Gaussianities in infiationary cosmology and a list of ref- 
erences). In the following we give a very brief summary 
of the analysis of non-Gaussianities in the matter bounce 
scenario. 

Non-Gaussianities are induced in any cosmological 
model simply because the Einstein equations are non- 
linear. In momentum space, the bispectrum contains 
amplitude and shape information. The bispectrum is a 
function of the three momenta. Momentum conserva- 
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tion implies that the three momenta have to add up to 
zero. However, this still leaves a rich shape information 
in the bispectrum in addition to the information about 
the overall amplitude. 

A formalism to compute the non-Gaussianities f or th e 
curvature fluctuation variable C was developed in jl29l |. 
Working in the interaction representation, the three- 
point function of ( is given to leading order by 



<C(i,fci)C(t,fc2)C(i,^3) > 



(120) 



= I dt' < [Cit,ki)C{t,k2)C{t,h),L„u{t')] > , 

where U corresponds to the initial time before which 
there are any non-Gaussianities. The square parenthe- 
ses indicate the commutator, and Lint is the interaction 
Lagrangian 

The interaction Lagrangian contains many terms. In 
particular, there are terms containing the time derivative 
of (. Each term leads to a particular shape of the bispec- 
trum. In an expanding universe such as in inflationary 
cosmology C = 0. However, in a contracting phase the 
time derivative of C does not vanish since the dominant 
mode is growing in time. Hence, there are new contri- 
butions to the shape which have a very different form 
from the shape of the terms which appear in inflationary 
cosmology. The larger value of the amplitude of the bis- 
pectrum follows again from the fact that there is a mode 
function which grows in time in the contracting phase. 

The three-point function can be expressed in the fol- 
lowing general form: 

<C(^i)C(^2)C(fc3)> = i27ry6iJ2h)^ 

xA{ki,k2,h) , (121) 

where ki = \ki\ and A is the shape function. In this 
expression we have factored out the dependence on the 
power spectrum . In inflationary cosmology it has be- 
come usual to express the bispectrum in terms of a non- 
Gaussianity parameter Jnl- However, this is only useful 
if the shape of the three point function is known. As a 
generalization, we here use |l27| 

IRI ft i: t \ lOMMjkM) f.ryry. 

\B\^^(ki,k2,k3) = — ^-p . (122) 

The computation of the bispectrum is tedious. In the 
case of the matter bounce (no entropy fluctuations) the 
result is 



-9j2kfk^ + 5j2k!k 



(123) 



-66 kfkpl+9 J2 kfkp^ 



This equation describes the shape which is predicted. 
Some of the terms (such as the last two) are the same 
as those which occur in single field slow-roll infiation, 
but the others are new. Note, in particular, that the new 
terms are not negligible. 

If we project the resulting shape function A onto some 
popular shape masks we 
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for the local shape (fci ^^2 = ^3). This is negative and 
of order 0(1). For the equilateral form (fci ~ k^ ~ k^) 
the result is 
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and for the folded form (fci — 2^2 = 2/03) one obtains the 
value 
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These amplitudes are close to what the Planck CMB 
satellite experiment will be able to detect. 

The matter bounce scenario also predicts a change 
in the slop e of the primordial power spectrum on small 
scales jl3Cll | : scales which exit the Hubble radius in the ra- 
diation phase retain a blue spectrum since the squeezing 
rate on scales larger than the Hubble radius is insuffi- 
cient to give longer wavelength modes a sufficient boost 
relative to the shorter wavelength ones. 



D. Problems of the Matter Bounce 

The main problem of the matter bounce scenario is 
that the physics which yields the non-singular bounce 
is not (yet) well established. It is clearly required to 
go beyond conventional quantum field theory coupled to 
General Relativity to obtain such a bounce. New Planck- 
scale physics needs to be invoked to obtain the required 
background cosmology. We remind the reader that it is 
precisely the absence of such new physics which needs 
to be invoked to argue for an inflationary background 
evolution. 

In terms of the trans-Planckian problem for cosmo- 
logical fluctuations, the bouncing cosmology has a clear 
advantage: the wavelength of the fluctuations which we 
are interested in always remains in the far infrared com- 
pared to the ultraviolet scales of the new physics. It can 
be shown that the correction terms due to the new ultra- 
violet physics on the evolution of fluctuations of interest 
to observers are negligible (see e.g. pslj). 

The second main problem of the matter bounce sce- 
nario may be the sensitivity of the bounce to assumptions 
on the initial conditions in the far past. The amplitude 
of the classical initial fluctuations must be small in or- 
der for the homogeneous background dynamics not to 
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be perturbed, and for the vacuum fluctuations to domi- 
nate on the scales relevant in the infrared (see e.g. (l3l| 
for a criticism of initial conditions in a related scenario, 
the Pre-Big-Bang scenario). In addition, the initial shear 
must be very small (in some bouncing cosmology back- 
grounds the bouncing solution is unstable to the smallest 
addition of shear. This, however, is not the case in other 
models such as the Horava-Lifshitz bounce). 

The matter bounce scenario does not address the cos- 
mological constant problem, but on the other hand the 
scenario does not depend on how the cosmological con- 
stant problem is solved. Note that inflationary cosmology 
is NOT robust in this respect. 



V. STRING GAS COSMOLOGY 

Another scenario of early universe cosmology which 
can provide an explanation for the current dat a is s tring 
gas cosmology. It is an "emergent universe" jl32| sce- 
nario in which the universe begins in a long hot and al- 
most static phase. The key input from string theory is 
that matter is a gas of closed fundamental strings com- 
pared to a gas of point particles as is assumed in Standard 
Big Bang cosmology. 



A. Principles and Background 

1. Principles 

String theory may be the best candidate we have at 
the present time for a quantum theory of gravity which 
unifies all four forces of nature. This theory is, however, 
currently not yet well understood beyond perturbation 
theory For applications to early universe cosmology, 
however, a non-perturbative understanding will be essen- 
tial. 

In the absence of a non-perturbative formulation of 
string theory, the approach to string cosmology which 
we have s uggested, str ing gas cosmology [l^l (see also 
[28| , and jlOl . Il34 Il35j | for reviews and more complete 
list of references), is to focus on symmetries and degrees 
of freedom which are new to string theory (compared to 
point particle theories) and which will be part of any non- 
perturbative string theory, and to use them to develop a 
new cosmology. The symmetry we make use of is T- 
duality, and the new degrees of freedom are the string 
oscillatory modes and the string winding modes. 

String gas cosmology is based on coupling a classical 
background which includes the graviton and the dilaton 
fields to a gas of closed strings (and possibly other ba- 
sic degrees of freedom of string theory such as "branes" 



Note that there are concrete proposals for a non- pert urbative 
formulation such as the AdS/CFT correspondence [133|| . 



[136j). All dimensions of space are taken to be compact, 
for reasons which will become clear later. For simplicity, 
we take all spatial directions to be toroidal and denote 
the radius of the torus by R. Strings have three types of 
states: momentum modes which represent the center of 
mass motion of the string, oscillatory modes which rep- 
resent the fluctuations of the strings, and winding modes 
counting the number of times a string wraps the torus. 

Since the number of string oscillatory states increases 
exponentially with energy, there is a limiting temperature 
for a gas of strings in thermal equilibrium, the Hagedorn 
temperature (29j Tjj- Thus, if we take a box of strings 
and adiabatically decrease the box size, the temperature 
will never diverge. This is the first indication that string 
theory has the potential to resolve the cosmological sin- 
gularity problem. 

The second key feature of string theory upon which 
string gas cosmology is based is T-duality. To introduce 
this symmetry, let us discuss the radius dependence of 
the energy of the basic string states: The energy of an 
oscillatory mode is independent of R, momentum mode 
energies are quantized in units of i.e. 

£;„ = nfij^ , (127) 

where Is is the string length and fi is the mass per unit 
length of a string. The winding mode energies are quan- 
tized in units of R, i.e. 

E,n = m^R, (128) 

where both n and m are integers. Thus, a new symme- 
try of the spectrum of string states emerges: Under the 
change 

R 1/R (129) 

in the radius of the torus (in units of Is) the energy spec- 
trum of string states is invariant if winding and momen- 
tum quantum numbers are interchanged 

(n, m) — ^ (m, n) . (130) 

The above symmetry is the simplest element of a larger 
symmetry group, the T-duality symmetry group which in 
general also mixes fluxes and geometry. The string vertex 
operators are consistent with this symmetry, and thus 
T-duality is a symmetry of perturbativc string theory. 
Postulating that T-duality extends to non-perturbative 
string theory leads [137J to the need of adding D-branes 
to the list of fundamental objects in string theory. With 
this addition, T-duality is expected to be a symmetry 
of non-perturbative string theory. Specifically, T-duality 
will take a spectrum of stable Type IIA branes and map it 
into a corresponding s pect rum of stable Type IIB branes 
with identical mass es Il38j |. 

As discussed in f27'|, the above T-duality symmetry 
leads to an equivalence between small a nd la,rge spaces, 
an equivalence elaborated on further in jl39l . ll4Cll | . 
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FIG. 12: The temperature (vertical axis) as a function of 
radius (iiorizontal axis) of a gas of closed strings in thermal 
equilibrium. Note the absence of a temperature singularity. 
The range of values of R for which the temperature is close to 
the Hagedorn temperature Th depends on the total entropy 
of the universe. The upper of the two curves corresponds to 
a universe with larger entropy. 
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1 > t 

p = p = rho/3 

FIG. 13: The dynamics of string gas cosmology. The vertical 
axis represents the scale factor of the universe, the horizon- 
tal axis is time. Along the horizontal axis, the approximate 
equation of state is also indicated. During the Hagedorn phase 
the pressure is negligible due to the cancellation between the 
positive pressure of the momentum modes and the negative 
pressure of the winding modes, after time the equation of 
state is that of a radiation-dominated universe. 



2. Background Cosmology 

That string gas cosmology will lead to a dynamical 
evolution of the early universe very different from what 
is obtained in standard and inflationary cosmology can 
already be seen by combining the basic ingredients from 
string theory discussed in the previous subsection. As 
the radius of a box of strings decreases from an initially 
very large value - maintaining thermal equilibrium - , the 
temperature first rises as in standard cosmology since the 
string states which are occupied (the momentum modes) 
get heavier. However, as the temperature approaches the 
Hagedorn temperature, the energy begins to flow into the 
oscillatory modes and the increase in temperature levels 
off. As the radius R decreases below the string scale, 
the temperature begins to decrease as the energy begins 
to flow into the winding modes whose energy decreases 
as R decreases (see Figure [T2]). Thus, as argued in p7j . 
the temperature singularity of early universe cosmology 
should be resolved in string gas cosmology. 

The equations that govern the background of string 
gas cosmology are not known. The Einstein equations 
are not the correct equations since they do not obey the 
T-duality symmetry of string theory. Many early studies 
of string gas cosm ology wer e based on using the dilaton 
gravity equations |l4ll4l43 |. However, these equations 



are not satisfactory, either. Firstly, we expect that string 
theory correction terms to the low energy effective ac- 
tion of string theory become dominant in the Hagedorn 
phase. Secondly, the dilaton gravity equations yield a 
rapidly changing dilaton during the Hagedorn phase (in 
the string frame). Once the dilaton becomes large, it be- 
comes inconsistent to focus on fundamental string states 
rather than brane states. In other words, using dilaton 



gravity as a background for string gas cosmology does not 
correctly reflect the S-duality symmetry of string theory. 
Recently, a background for string g as co smology includ- 
ing a rolling tachyon was proposed jl44l | which allows a 
background in the Hagedorn phase with constant scale 
factor and consta nt di laton. Another study of this prob- 
lem was given in |l45| . 

Some conclusions about the time-temperature relation 
in string gas cosmology can be derived based on thermo- 
dynamical considerations alone. One possibility is that 
R starts out much smaller than the self-dual value and 
increases monotonically. From Figure [T^] it then follows 
that the time-temperature curve will correspond to that 
of a bouncing cosmology. Alternatively, it is possible that 
the universe starts out in a meta-stable state near the 
Hagedorn temperature, the Hagedorn phase, and then 
smoothly evolves into an expanding phase dominated by 
radiation like in standard cosmology (Figure fTS]). Note 
that we are assuming that not only the scale factor but 
also the dilaton is constant in time. 

The transition between the quasi-static Hagedorn 
phase and the radiation phase at the time tn is a conse- 
quence of the annihilation of string winding modes into 
string loops (see Figure [T4| . Since this process corre- 
sponds to the production of radiation, we denote this 
time by the same symbol t/j as the time of reheating in 
inflationary cosmology. As pointed out in [23] ; this an- 
nihilation process only is possible in at most three large 
spatial dimensions. This is a simple dimension counting 
argument: string world sheets have measure zero inter- 
section probability in more than four large space-time 
dimensions. Hence, string gas cosmology may provide a 
natural mechanism for explaining why there are exactly 
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FIG. 14: The process by which string loops are produced via 
the intersection of winding strings. The top and bottom lines 
are identified and the space between these lines represents 
space with one toroidal dimension un-wrapped. 



three large spatial dimensions. This argument was sup- 
ported by numerical studies of st ring evolut ion i n three 
and four spatial dimensions |146{ (see also jl47i |). The 
flow of energy from winding modes to string l oops can be 
modelled by effective Boltzmann equations | 148| analo- 
gous to those used to describe the flow of energy between 
infinite cosmic strings and cosmic string loops (see e.g. 
[ii-ls^l for reviews). 

Several comments are in place c oncer ning the above 
mechanism. First, in the analysis of jl48j | it was assumed 
that the string interaction rates were time-independent. 
If the dynamics of the Hagedorn phase is modelled by 
dilaton gravity, the dilaton is rapidly changing during 
the phase in whi ch the stri ng frame scale factor is static. 
As discussed in [l49l Il50l | (see also (l5l( ). in this case 
the mechanism which tells us that exactly three spatial 
dimensions become macroscopic does not work. Another 
comment concerns the isotropy of the three large dimen- 
sions. In contrast to the situation in Standard cosmol- 
ogy, in string gas cosm ology the anisotropy decreases 
in the expanding phase [152| . Thus, there is a natural 
isotropization mechanism for the three large spatial di- 
mensions. 



3. Moduli StaHUzation 

In the following, we shall assume that either the mech- 
anism of [27] for setting in motion the preferential expan- 
sion of exactly three spatial dimensions works, or, alter- 
natively, that three dimensions are distinguished from the 
beginning as being large. In either case, we must address 
the moduli stabilization problem, i.e. we must show that 
the radii (radions) and shapes of the extra dimensions are 
stabilized. This is a major challenge in string-motivated 
field theory models of higher dimensions. The situation 
in string gas cosmology is much better in comparison: all 
moduli except for the dilaton are stabilized without the 
need of introducing extra ingredients such as fluxes or 
special non-perturbative effects. 

Radion stab ilization in the string frame was initially 
studied in |l53j | . The basic idea is that the winding modes 
about the extra spatial dimensions provide a confining 
force which prevents the radii from increasing whereas 



the momentum modes provide a force which resists the 
complete contraction. Thus, there will be a stable mini- 
mum of the effective potential for the radion. 

In order to make contact with late time cosmology, it 
is important to consider the issue of radion stabilization 
when the dilaton is frozen, or, mor e ge i ieral ly, in the 
Einstein frame. As was discussed in jl54 Il55| | (see also 
earlier comments in [l53j |). the existence of string modes 
which are massless at the self-dual radius is crucial in 
obtaining radion stabilization in the Einstein frame (for 
more general studies of t he impor tance of massless modes 
in string cosmology see jl56lll57l |). Such massless modes 
do not exist in all known string theories. They ex ist in 
the Heterotic theory, but not in Type II th eories |l37l | . 
The following discussion is taken from |l58j |. 

Let us consider the equations of motion which arise 
from coupling the Einstein action to a string gas. In the 
anisotropic setting when the metric is taken to be 

6 

ds^ = d<2 _ a{tfdy? - baitfdyl , (131) 

where the i/a are the internal coordinates, the equation 
of motion for ba becomes 

b^ + (3H+ y ^)b^ ^ y 87tG-(^^S (132) 

where ^m,n is the number density of (m, n) strings, e^.n 
is the energy of an individual (m, n) string, and g is the 
determinant of the metric. The source term S depends on 
the quantum numbers of the string gas, and the sum runs 
over all momentum and winding number vectors m and 
n, respectively (note that n and m are six-vectors, one 
component for each internal dimension). If the number 
of right- moving oscillator modes is given by N , then the 
source term for fixed m and n vectors is 

S = y{'^f^ynlbl+^^[{n,n) + {n.m)+2{N^l 

(133) 

To obtain this equation, we have made use of the mass 
spectrum of string states and of the level matching con- 
ditions. In the case of the bosonic superstring, the mass 
M of a string state with fixed m,n,N and N, where N 
and N are the number of right- and left-moving oscillator 
states, on a six-dimensional torus whose radii are given 
by ba is 

= {^f-y nlbl + 2(iV + TV - 2) , (134) 

a 

and the level matching condition reads 

N = {n,m) + N, (135) 

where (n, m) indicates the scalar product of n and m in 
the trivial metric. 



29 



There are modes which are massless at the self-dual 
radius ha = 1. One such mode is the graviton with n = 
TO = and = 1. The modes of interest to us are modes 
which contain winding and momentum, namely 

• = 1, (to, to) = 1, {m,n) — —1 and {n,n) — 1; 

• A'^ = 0, (to, to) = 1, (to, n) = 1 and (n, n) — 1; 

• N = (to, m) = 2, (m, n) = and (n, n) — 2. 

Note that these modes survive in t he H eterotic string 
theory, but do not survive the GSO |137| truncation in 
Type II string theories. 

In string theories which admit massless states (i.e. 
states which are massless at the self-dual radius), these 
states will dominate the initial partition function. The 
background dynamics will then also be dominated by 
these states. To understand the effect of these strings, 
consider the equation of motion (|132[) with the source 
term (|133p . The first two terms in the source term cor- 
respond to an effective potential with a stable minimum 
at the self-dual radius. However, if the third term in 
the source does not vanish at the self-dual radius, it will 
lead to a positive potential which causes the radion to 
increase. Thus, a condition for the stabilization of ba at 
the self-dual radius is that the third term in (|133p van- 
ishes at the self-dual radius. This is the case if and only 
if the string state is a massless mode. 

The massless mode s hav e other nice features which are 
explored in detail in |l55|. They act as radiation from 
the point of view of our three large dimensions and hence 
do not lead to a over-abundance problem. As our three 
spatial dimensions grow, the potential which confines the 
radion becomes shallower. However, rather surprisingly, 
it turns out the the potential remains steep enough to 
avoid fifth force constraints. 

In the presence of massless string states, the shape 
moduli also can be stabilized, at leas t in the simple 
toroidal backgrounds considered so far [l59j . To study 
this issue, we consider a metric of the form 



- Gmndy^'dy^ , 



(136) 



where the metric of the internal space (here for simplic- 
ity considered to be a two-dimensional torus) contains a 
shape modulus, the angle between the two cycles of the 
torus: 



Gil — G 
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and 



Gi2 = G21 = sine, 



(137) 



(138) 



where 6* = corresponds to a rectangular torus. The 
ratio between the two toroidal radii is a second shape 
modulus. However, from the discussion of the previous 
subsection we already know that each radion individu- 
ally is stabilized at the self-dual radius. Thus, the shape 
modulus corresponding to the ratio of the toroidal radii 



is fixed, and the angle is the only shape modulus which 
is yet to be considered. 

Combining the 00 and the 12 Einstein equations, we 
obtain a harmonic oscillator equation for 9 with = as 
the stable fixed point. 



(139) 



where K is a, constant whose value depends on the quan- 
tum numbers of the string gas. In the case of an expand- 
ing three-dimensional space we would have obtained an 
additional damping term in the above equation of mo- 
tion. We thus conclude that the shape modulus is dy- 
namically stabilized at a value which maximizes the area 
to circumference ratio. 

The only modulus which is not stabilized by string 
winding and momentum modes is the dilaton. Recently, 
it has been show 160 ] that a gaugino condensation mech- 
anism (similar to those used in string inflation model 
building) can be introduced which generates a stabiliz- 
ing potential for the dilaton without interfering with the 
radion stabilization force provided by the string winding 
and momentum modes. 

In the next subsection we turn to the predictions of 
string gas cosmology for observations. These predictions 
do not depend on the details of the theory, but only on 
three inputs. The first is the existence of a quasi-static 
initial phase in thermal equilibrium. The second condi- 
tion is the applicability of the Einstein field equations 
for fluctuations on infrared scales (scales of the order of 
1mm), many orders of magnitude larger than the string 
scale, and the third is a holographic scaling of the specific 
heat capacity. The role of the last condition will become 
manifest below. 



B. Fluctuations in String Gas Cosmology 

1. Overview 

At the outset of this section, let us recall the mecha- 
nism by which inflationary cosmology leads to the pos- 
sibility of a causal generation mechanism for cosmolog- 
ical fluctuations which yields an almost scale-invariant 
spectrum of perturbations. The space-time diagram of 
inflationary cosmology is sketched in Figure [S] 

During the period of inflation, the Hubble radius 



lH{t) 



a 
a 



(140) 



is approximately constant. In contrast, the physical 
length of a fixed co-moving scale (labelled by k in the fig- 
ure) is expanding exponentially. In this way, in inflation- 
ary cosmology scales which have microscopic sub-Hubble 
wavelengths at the beginning of inflation are red-shifted 
to become super-Hubble-scale fluctuations at the end of 
the period of inflation. After inflation, the Hubble radius 
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increases linearly in time, faster than the physical wave- 
length corresponding to a fixed co-moving scale. Thus, 
scales re-enter the Hubble radius at late times. 

Since inflation red-shifts any classical fluctuations 
which might have been present at the beginning of the 
inflationary phase, fluctuations in inflationary cosmology 
are generated by quantum vacuum perturbations. The 
fluctuations begin in their quantum vacuum state at the 
onset of inflation. Once the wavelength exceeds the Hub- 
ble radius, squeezing of the wave-function of the fluc- 
tuations sets in (see [S^ HI])- This squeezing plus the 
de-coherence of the fluctuations due to the interaction 
between short and long wavelength modes generated by 
the intrinsic non-linearities in both the gravitational and 
matter sectors of the theory (see i70] for recent dis- 
cussions of this aspect and references to previous work) 
lead to the classicalization of the fluctuations on super- 
Hubble scales. 

Let us now turn to the cosmological background of 
string gas cosmology represented in Figure 1131 This 
string gas cosmology background yields the space-time 
diagram sketched in Figure [131 As in Figure [5l the verti- 
cal axis is time and the horizontal axis denotes the phys- 
ical distance. For times t < t^, we are in the static 
Hagedorn phase and the Hubble radius is infinite. For 
t > tii, the Einstein frame Hubble radius is expanding as 
in standard cosmology. The time is when the string 
winding modes begin to decay into string loops, and the 
scale factor starts to increase, leading to the transition 
to the radiation phase of standard cosmology. 

Let us now compare the evolution of the physical wave- 
length corresponding to a fixed co-moving scale with that 
of the Einstein frame Hubble radius H^^{t). The evolu- 
tion of scales in string gas cosmology is identical to the 
evolution in standard and in inflationary cosmology for 
t > tji. If we follow the physical wavelength of the co- 
moving scale which corresponds to the current Hubble 
radius back to the time tij, then - taking the Hagedorn 
temperature to be of the order 10^^ GeV - we obtain a 
length of about 1 mm. Compared to the string scale and 
the Planck scale, this is in the far infrared. 

The physical wavelength is constant in the Hagedorn 
phase since space is static. But, as we enter the Hagedorn 
phase going back in time, the Hubble radius diverges to 
inflnity. Hence, as in the case of inflationary cosmology, 
fluctuation modes begin sub- Hubble during the Hagedorn 
phase, and thus a causal generation mechanism for fiuc- 
tuations is possible. 

However, the physics of the generation mechanism is 
very different. In the case of inflationary cosmology, 
fluctuations are assumed to start as quantum vacuum 
perturbations because classical inhomogeneities are red- 
shifting. In contrast, in the Hagedorn phase of string 
gas cosmology there is no red-shifting of classical matter. 
Hence, it is the fluctuations in the classical matter which 
dominate. Since classical matter is a string gas, the domi- 
nant fluctuations are string thermodynamic fiuctuations. 

Our proposal for string gas structure formation is the 
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FIG. 15: Space-time diagram (sketch) showing the evolution 
of fixed co-moving scales in string gas cosmology. The vertical 
axis is time, the horizontal axis is physical distance. The solid 
curve represents the Einstein frame Hubble radius which 
shrinks abruptly to a micro-physical scale at Ir and then in- 
creases linearly in time for t > ta. Fixed co-moving scales 
(the dotted lines labeled by fci and ^2) which are currently 
probed in cosmological observations have wavelengths which 
are smaller than the Hubble radius before tn. They exit the 
Hubble radius at times ti{k) just prior to tn, and propagate 
with a wavelength larger than the Hubble radius until they 
reenter the Hubble radius at times tf{k). 



following [sOl (see [l6l j for a more detailed descrip- 
tion). For a fixed co- moving scale with wavenumber k 
we compute the matter fluctuations while the scale in 
sub-Hubble (and therefore gravitational effects are sub- 
dominant). When the scale exits the Hubble radius at 
time ti(k) we use the gravitational constraint equations 
to determine the induced metric fluctuations, which are 
then propagated to late times using the usual equations 
of gravitational perturbation theory. Since the scales we 
are interested in are in the far infrared, we use the Ein- 
stein constraint equations for fluctuations. 



31 



2. Spectrum of Cosmological Fluctuations 

We write the metric including cosmological perturba- 
tions (scalar metric fluctuations) and gravitational waves 
in the following form (using conformal time rf) 

ds^ = a^(??) {(1 + 2^)dT]^ - [(1 - 2$)% + h,j]dx'dx^} . 

(141) 

We have fixed the gauge (i.e. coordinate) freedom for the 
scalar metric perturbations by adopting the longitudinal 
gauge and we have taken matter to be free of anisotropic 
stress. The spatial tensor /iy (x,t) is transverse and 
traceless and represents the gravitational waves. 

Note that in contrast to the case of slow-roll infla- 
tion, scalar metric fluctuations and gravitational waves 
are generated by matter at the same order in cosmo- 
logical perturbation theory. This could lead to the ex- 
pectation that the amplitude of gravitational waves in 
string gas cosmology should be generically larger than 
in inflationary cosmology. This expectation, however, is 
not realized 162] since there is a different mechanism 
which suppresses the gravitational waves relative to the 
density perturbations (namely the fact that the gravi- 
tational wave amplitude is set by the amplitude of the 
pressure, and the pressure is suppressed relative to the 
energy density in the Hagedorn phase). 

Assuming that the fluctuations are described by the 
perturbed E instein eq uations (they are not if the dilaton 
is not fixed (l63l llG^ V then the spectra of cosmological 
perturbations 3" and gravitational waves h are given by 
the e nergy-momentum fluctuations in the following way 

{mk)\^) - 16n^ G^k-^{ST%{k)dT%ik)) , (142) 

where the pointed brackets indicate expectation values, 
and 

(|/i(fc)|2) ^ 16n'^ G'^k-'^{ST' J {k)5T' J (k)) , (143) 



where on the right hand side of (|143p we mean the average 
over the correlation functions with i ^ j, and h is the 
amplitude of the gravitational waves . 

Let us now use p42p to determine the spectrum of 
scalar metric fluctuations. We first calculate the root 
mean square energy density fluctuations in a sphere of 
radius R = k~^. For a system in thermal equilibrium 
they are given by the specific heat capacity Cy via 



of c losed strings. This computation was carried out 



gas 

in Il65l with the result 



Cv 



T (1 - T/Th) 



(145) 



The specific heat capacity scales holographically with the 
size of the box. This result follows rigorously from evalu- 
ating the string partition function in the Hagedorn phase. 
The result, however, can also be understood heuristically: 
in the Hagedorn phase the string winding modes are cru- 
cial. These modes look like point particles in one less 
spatial dimension. Hence, we expect the specific heat ca- 
pacity to scale like in the case of point particles in one 
less dimension of space . 

With these results, the power spectrum P{k) of scalar 
metric fiuctuations can be evaluated as follows 



-P*(fc) 



1 

2^ 



■A:3|$(fc)|2 



(146) 



= 8G^fc-i < |(5/9(fc)p > . 
= SG^fc^ < {5Mf >R 
= 8G^k-^ < {6pf >R 
T 1 



= 8G 



1 - T/Th 



where in the first step we have used p42p to replace the 
expectation value of |<I>(A:)p in terms of the correlation 
function of the energy density, and in the second step we 
have made the transition to position space 

The first conclusion from the result (jl46|) is that the 
spectrum is approximately scale-invariant {P{k) is inde- 
pendent of fc). It is the 'holographic' scaling Cv{R) ~ 
which is responsible for the overall scale-invariance of the 
spectrum of cosmological perturbations. However, there 
is a small k dependence which comes from the fact that 
in the above equation for a scale k the temperature T 
is to be evaluated at the time ti{k). Thus, the factor 
(1 — T/Th) in the denominator is responsible for giving 
the spectrum a slight dependence on k. Since the tem- 
perature slightly decreases as time increases at the end 
of the Hagedorn phase, shorter wavelengths for which 
ti{k) occurs later obtain a smaller amplitude. Thus, the 
spectrum has a slight red tilt. 



3. Key Prediction of String Gas Cosmology 



(144) 



The specific heat of a gas of closed strings on a torus of 
radius R can be derived from the partition function of a 



As discovered in jl62|, the spectrum of gravitational 
waves is also nearly scale invariant. However, in the 
expression for the spectrum of gravitational waves the 



The gravitational wave tensor hij can be written as the ampli- 
tude h multiplied by a constant polarization tensor. 



We emphasize that it was important for us to have compact 
spatial dimensions in order to obtain the winding modes which 
are crucial to get the holographic scaling of the thermodynamic 
quantities. 
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factor (1 — T/Th) appears in the numerator, thus lead- 
ing to a shght blue tilt in the spectrum. This is a pre- 
diction with which the cosmological effects of string gas 
cosmology can be distinguished from those of inflation- 
ary cosmology, where quite generically a slight red tilt for 
gravitational waves results. The physical reason is that 
large scales exit the Hubble radius earlier when the pres- 
sure and hence also the off-diagonal spatial components 
of Tf^^ are closer to zero. 

Let us analyze this issue in a bit more detail and 
estimate the dimensionless power spectrum of gravita- 
tional waves. First, we make some general comments. In 
slow-roll inflation, to leading order in perturbation the- 
ory matter fluctuations do not couple to tensor modes. 
This is due to the fact that the matter background field 
is slowly evolving in time and the leading order gravita- 
tional fluctuations are linear in the matter fluctuations. 
In our case, the background is not evolving (at least at 
the level of our computations), and hence the dominant 
metric fluctuations are quadratic in the matter field fluc- 
tuations. At this level, matter fluctuations induce both 
scalar and tensor metric fluctuations. Based on this con- 
sideration we might expect that in our string gas cosmol- 
ogy scenario, the ratio of tensor to scalar metric fluctu- 
ations will be larger than in simple slow-roll inflationary 
models. However, since the amplitude h of the gravi- 
tational waves is proportional to the pressure, and the 
pressure is suppressed in the Hagedorn phase, the am- 
plitude of the gravitational waves will also be small in 
string gas cosmology. 

The method for calculating the spectrum of gravita- 
tional waves is similar to the procedure outlined in the 
last section for scalar metric fluctuations. For a mode 
with fixed co-moving wavenumber k, we compute the cor- 
relation function of the off-diagonal spatial elements of 
the string gas energy- momentum tensor at the time U (k) 
when the mode exits the Hubble radius and use (|143|) 
to infer the amplitude of the power spectrum of gravita- 
tional waves at that time. We then follow the evolution of 
the gravitational wave power spectrum on super-Hubble 
scales for t > ti{k) using the equations of general rela- 
tivistic perturbation theory. 

The power spectrum of the tensor modes is given by 
(fT43)) : 

Phik) ^ 16t:^ G^k-'^k^ {ST' J {k)ST' J (k)) (147) 

for i ^ j. Note that the k^ factor multiplying the mo- 
mentum space correlation function of T'j gives the po- 
sition space correlation function C j'^ j{R) , namely the 
root mean square fluctuation of in a region of radius 
R — k~^ (the reader who is skeptical about this point is 
invited to check that the dimensions work out properly). 
Thus, 



Ph{k) = 167T^G^k^''C'j'j{R) 



(148) 



The correlation function C*/j on the right hand side of 
the above equation follows from the thermal closed string 



partition function and was computed explicitly in \l6t 
(see also jT3|)- We obtain 



Ph{k) ^ m7T^G^^{l-T/TH)\n 



(149) 

which, for temperatures close to the Hagedorn value re- 
duces to 



Ph{k) 
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P^k' 



■{1-T/Th] 



(150) 

This shows that the spectrum of tensor modes is - to 
a first approximation, namely neglecting the logarith- 
mic factor and neglecting the fc-dependence of T{ti{k)) - 
scale-invariant. 

On super-Hubble scales, the amplitude h of the grav- 
itational waves is constant. The wave oscillations freeze 
out when the wavelength of the mode crosses the Hub- 
ble radius. As in the case of scalar metric fluctuations, 
the waves are squeezed. Whereas the wave amplitude re- 
mains constant, its time derivative decreases. Another 
way to see this squeezing is to change variables to 



V'(?;,x) — a(ri)h(ri,x.) 



(151) 



in terms of which the action has canonical kinetic term. 
The action in terms of ip becomes 



(152) 



from which it immediately follows that on super-Hubble 
scales ij) ^ a. This is the squeezing of gravitational waves 

Since there is no /c-dependence in the squeezing fac- 
tor, the scale-invariance of the spectrum at the end of 
the Hagedorn phase will lead to a scale-invariance of the 
spectrum at late times. 

Note that in the case of string gas cosmology, the 
squeezing factor 2(77) does not differ substantially from 
the squeezing factor a{rii) for gravitational waves. In 
the case of inflationary cosmology, 2(77) and 0(77) dif- 
fer greatly during reheating, leading to a much larger 
squeezing for scalar metric fluctuations, and hence to a 
suppressed tensor to scalar ratio of fluctuations. In the 
case of string gas cosmology, there is no difference in 
squeezing between the scalar and the tensor modes. 

Let us return to the discussion of the spectrum of grav- 
itational waves. The result for the power spectrum is 
given in (|150p . and we mentioned that to a first approx- 
imation this co rresp onds to a scale-invariant spectrum. 
As realized in |l62| . the logarithmic term and the k- 
dependence of T(ti{k)) both lead to a small blue-tilt of 
the spectrum. This feature is characteristic of our sce- 
nario and cannot be reproduced in inflationary models. 
In inflationary models, the amplitude of the gravitational 
waves is set by the Hubble constant H. The Hubble con- 
stant cannot increase during inflation, and hence no blue 
tilt of the gravitational wave spectrum is possible. 
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A heuristic way of understanding the origin of the 
shght blue tih in the spectrum of tensor modes is as fol- 
lows. The closer we get to the Hagedorn temperature, 
the more the thermal bath is dominated by long string 
states, and thus the smaller the pressure will be com- 
pared to the pressure of a pure radiation bath. Since the 
pressure terms (strictly speaking the anisotropic pressure 
terms) in the energy-momentum tensor are responsible 
for the tensor modes, we conclude that the smaller the 
value of the wavenumber k (and thus the higher the tem- 
perature T{ti{k)) when the mode exits the Hubble radius, 
the lower the amplitude of the tensor modes. In contrast, 
the scalar modes are determined by the energy density, 
which increases at T{ti{k)) as k decreases, leading to a 
slight red tilt. 

The reader may ask about the predictions of string 
gas cosmology for non-Gaussianities. The answer is 
|167l | that the non-Gaussianities from the thermal string 
gas perturbations are Poisson-suppressed on scales larger 
than the thermal wavelength in the Hagedorn phase. 
However, if the spatial sections are initially larg e, th en it 
is possible that a network of cosmic superstrings [168| will 
be left behind. These strings - if stable - would achieve 
a scaling solution (constant number of strings crossing 
each Hubble volume at each time [85l-l87l|). Such strings 
give rise to linear discontinuities in the CMB tempera- 
ture maps |l69l |. lines which can be searched for using 
edge detection algorithms such as the Canny algorithm 
(sec [170.] for recent feasibility studies). 



C. Problems of String Gas Cosmology 

The main problem of the current implementation of 
string gas cosmology is that it does not provide a quan- 
titative model for the Hagedorn pha se. The criticisms 

of string gas cosmology raised in |l64l . Il7l| are based on 
assuming that dilaton gravity should be the background. 
However, as mentioned earlier, this is not a reasonable as- 
sumption since the coupling of string gas matter to dila- 
ton gravity is not S-duality invariant. The quasi-static 
early Hagedorn phase should also have constant dilaton. 
In addition, we do not expect that at high densities such 
as the Hagedorn density an simple effective action such as 
that of dilaton gravity will apply. Nevertheless, to make 
the string gas cosmology scenario into a real theory, it 
is crucial to obtain a good understanding of the back- 
gro und d ynamics. For some initial steps in this direction 
see pi] . 

A second problem of string gas cosmology is the size 
problem (and the related entropy problem). If the string 
scale is about lO^^GeV as is preferred in early heterotic 
superstring models, then the radius of the universe during 
the Hagedorn phase must be many orders of magnitude 
larger than the string scale. Without embedding string 
gas cosmology into a bouncing cosmology, it seems un- 
natural to demand such a large initial size. This problem 
disappears if the Hagedorn phase is preceded by a phase 



of contraction, as in the model of [172| . In this case, 
however, it is non-trivial to arrange that the Hagedorn 
phase lasts sufhciently long to maintain thermal equilib- 
rium over the required range of scales. 

It should be noted, however, that some of the con- 
ceptual problems of inflationary cosmology such as the 
trans-Planckian problem for fluctuations, do not arise in 
string gas cosmology. As in the case of the matter bounce 
scenario, the basic mechanism of the scenario is insensi- 
tive to what sets the cosmological constant to its observed 
very small value. 

VI. CONCLUSIONS 

The main messages of these lectures are the following: 

• It is possible to explore physics of the very early 
universe using current cosmological observations. 
Given the large amount of new data which is ex- 
pected over the next decade, early universe cosmol- 
ogy will be a vibrant held of research. The informa- 
tion about the very early universe is transferred to 
the current time mainly via imprints on the spec- 
trum of cosmological fluctuations. 

• The theory of cosmological fluctuations is the key 
tool of modern cosmology. It links current obser- 
vations with early universe cosmology. The theory 
can be applied to any background cosmology, not 
just inflationary cosmology. 

• The Hubble radius plays a key role in the evo- 
lution of fluctuations. On sub-Hubble scales the 
microphysical forces dominate, whereas on super- 
Hubble scales matter forces freeze out and gravity 
dominates. In order to have a causal mechanism 
of structure formation, it is therefore crucial that 
scales of current interest in cosmology originate at 
very early times inside the Hubble radius, and that 
they then propagate over an extended period of 
time on super-Hubble scales. 

• Inflationary cosmology provides a mechanism to ex- 
plain the origin of structure in the universe via 
causal physics, but it is not the only one. In these 
lectures we have presented two alternative scenar- 
ios, the "matter bounce" and "string gas cosmol- 
ogy" • 

• Inflationary cosmology was predictive: it predicted 
the detailed shape of the angular power spectrum 
of CMB anisotropies - more than 15 years before 
the observations. The challenge for alternative sce- 
narios is to identify clean predictions with which 
the new scenarios can be distinguished from those 
of inflation. In the case of the matter bounce the 
prediction we have identified is a particular shape 
of the bispectrum, in the case of string gas cosmol- 
ogy the key prediction is a slight blue tilt in the 
spectrum of gravitational waves. 



34 



Inflationary cosmology has been under development for 
30 years. The framework is based on Einstein gravity 
coupled to classical scalar field matter and is well es- 
tablished. However, when applied to early universe cos- 
mology there are important conceptual problems which 
arise. In contrast, the "matter bounce" and "string gas" 
scenarios of structure formation are much more recent. 
They both must involve physics which goes beyond clas- 
sical field theory coupled to Einstein gravity. To come 
up with a consistent model is a major challenge which 
has not yet been satisfactorily met. There is a lot of 
interesting work ahead of us. 
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